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We consider the simplest gauge theories given by one- and two- matrix integrals and concentrate on 
. their stringy and geometric properties. We remind general integrable structure behind the matrix 

integrals and turn to the geometric properties of planar matrix models, demonstrating that they are 
universally described in terms of integrable systems directly related to the theory of complex curves. 
We study the main ingredients of this geometric picture, suggesting that it can be generalized beyond 
one complex dimension, and formulate them in terms of the quasiclassical integrable systems, solved 
by construction of tau-functions or prepotentials. The complex curves and tau-functions of one- and 
two- matrix models are discussed in detail. 
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1 Gauge fields, geometry and string theory 

The relation between gauge fields and strings is still one of the most interesting problems of modern theoretical 
physics, see e.g. £Q. At present it goes far beyond the well-known fact, that spectrum of open strings in flat back- 
ground contains the vector field, which couples to the simplest open string excitation, and reparameterization 
invariance of the quantum mechanics of string modes requires it to be massless. 

The study of moduli space of closed string theories suggests that even there it is natural to introduce open 
strings, ending only on certain hypersurfaces or D-branes 2'. The stack of N D-branes leads in this way to 
appearence of the matrix or non-Abelian gauge field, where the indices of the matrix label the particular branes 
from the stack so that the corresponding string starts from and ends on. This is a geometric way to formulate 
the SU(N) Yang-Mills theory, and one finds in this way that the gauge theory, or the stack of D-branes, modifies 
the background geometry of the closed string theory, see detailed discussion of these issues, say in [HUH EI]. This 
allows to conjecture that gauge theory can be reformulated in terms of certain gravity picture, presumably in 
some multi-dimensional target space. 

The non-Abelian gauge theory must not have colored states (confinement). In close string theory language 
it means that, when closed-string background is modified, the open sting states should disappear from the 
spectrum. In terms of partition functions (or generating functions for the string amplitudes) this requirement 
acquires the form of equality 

-Zgaugc = exp (Faring ) (1.1) 

* Based on lectures presented at several schools on mathematical physics and the talks at "Complex geometry and string theory" 
and the Polivanov memorial seminar. 
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where Z gauge is statistical sum in some gauge theory, while its free energy String should be a partition function 
of a certain closed string theory. 

The most simple (but not the only!) geometric interpretation of the relation is the -^-expansion |S] 

in the inverse size of matrices from the side of matrix gauge theory. The perturbation series in quantum field 
theory has vanishing radius of convergency, if one sums over all possible diagrams. However, in matrix theory 
the Feynman graphs are fat, and each Feynman diagram can be uniquely assigned a two-dimensional surface of 
particular topology where it can be drawn without self intersections. Hence, an extra parameter - the genus of 
graph or of such two dimensional Riemann surface appears - and the series of the graphs of given topology can 
be possibly summed up within some finite radius of convergency. 

This means that for any matrix theory of the N x N matrices one gets in fact a " double" expansion 

oo 

String =J2 N2 ~ 29 E F 9> hXk (!- 2 ) 
3=0 h 

where g is the introduced genus of Riemann surface, corresponding the Feynman graph with V vertices E 
edges and F faces with the Euler characteristics \g — V — E + F = 2 — 2g. An extra sum over the holes 
h = 2g — 2 + F or faces of graphs, drawn on a given Riemann surface, can be thought of as summing (with a 
generating parameter X to be identified slightly later) over all contributions of open the strings, which, being 
summed up, disappear from the spectrum of the theory, modifying the geometry og closed strings. 

A sample example of gauge theory we study below is the zero-dimensional theory, or the matrix integral 
(considered, first, in this context in [7]) 

Z = J dfce-i™^*) (1.3) 
where $ is N x N matrix (for example, Hermitian), and W($) is some generic polynomial 

W(<S>)=J2 tk ® k (1.4) 

fc>0 

For the matrices of finite size N x N this partition function is just 7V 2 -tuple integral, which nevertheless 
possesses almost all known nontrivial properties of the quantum gauge theory path integrals. Obviously the 
"action" TrW($) in l|1.3fl is invariant under the zero-dimensional gauge transformations (where the derivative 
term is absent) $ — > W&U. The integration measure in (|1.3fl is just a certainly normalized product of the 
differentials of all matrix elements 

where, as in any gauge theory, the normalization factor Vn contains the volume of " orbit" of the unitary group 
U(N)/U(1) N together with the "eigenvalue-permutation" factor TV!. This normalization becomes natural after 
separation of the physical and gauge degrees of freedom - the diagonalization of the matrix $ by unitary gauge 
transformation $ = U^tftU, where <j> is diagonal. The volume of the orbit U{N)/U(\) N is then cancelled from 
denominator after trivial integration over the gauge degrees of freedom. 

Taking, as an example the potential W(3>) = + £^-$ 3 + . . . (where in analogy with the Yang-Mills 
theory we denoted the coupling constant by gyAi), introducing the 't Hooft coupling X — g\ M N and expanding 
the integral (|1 . 3|) after rescaling of the field one gets exactly expansion of the form 1)1. 2[l . since for a graph 
with V vertices, E propagators and H closed loops or holes the corresponding contribution is weighted by 

E\ V f±) E N h = n v-e+h x e-v (L6) 
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This counting is universal and does not depend on particular model, i.e. the same is true for realistic four- 
dimensional gauge theory with an SU(N) gauge group (direct generalization of phenomenological SU(2) or 
SU{3)). 

Now, as we already mentioned, the series in X over the contribution of open strings (at fixed genus g\) can 
be in principle summed up. This leads to the formula 

oo 

String =]T^ 2 - 2S F 9 (A) 

g=o 

where F g can be thought of as genus-g partition functions of some closed string theory. More strictly, it means 
that parameter X can be interpreted as a parameter of closed string background. It is very important to point 
out that string interpretation requires smooth Riemann surfaces in (|1.7|l . achieved only in the regime when 
contribution of all gauge diagrams is essential, i.e. when X 3> 1 or in the non-perturbative regime for the gauge 
theory. It is also necessary to point out that there is no direct geometric language, describing in this way gauge 
theory at weak coupling, since the corresponding limit on string theory side is very singular 1 . 

This is both good and bad. Good, since the gauge/string duality conjecture allows in principle, if we trust 
it, to predict non-trivial effects in the non-perturbative phase of gauge theory, i.e. in the regime which is mostly 
interesting for present theory of elementary particles. And it is bad, since it becomes very hard to test the 
conjecture itself, since it is not easy to provide good examples, for which both sides of equality can be 

calculated by independent methods. Basically it is the reason, why study of the matrix integrals (|1 is so 
important - since this is the simplest test example of the gauge/string correspondence. As we see below, already 
in this case it is a nontrivial problem, with lots of nice structures hidden inside, which are not only interesting by 
themselves, but, as we try to demonstrate in the very end, can be also applicable to testing higher-dimensional 
examples of the gauge/string correspondence. 

The -^--expansion can be also understood as quasiclassical computation of the matrix integral (|1.3fl . where 
the square of gauge coupling g\ M — h can be identified with the closed string coupling ti = g s or the Planck 
constant, which governs an expansion over the closed string loops. In what follows it is useful to rename the 
't Hooft coupling X, introducing the variable 

t = HN (1.8) 

which, in particular, will remain fixed and finite under quasiclassical limit N — > oo together with h — > 0. The 
■^-expansion l|1.7|l then can be equivalently written as a quasiclassical one 

oo 

String = ^T d h 2a - 2 F g {t ) (1.9) 
9=0 

where the coefficients differ from (f 1 . T|l by different scaling in to = X. 

Example: the gaussian case 

For the potential W / ( < i>) = ^Q 2 — ti<& the integral l|1.3fl can be immediately computed: 

Zgauss = J d$ e -* Tr (^ 2 -^) e -4# J dc&e-^* 2 = C N e~^ (2^)* (1.10) 

with Cn — y^- The A-dependence of the volume of unitary group ensure for (|1.1U|) that both J^r log Z gauss (N) 
and the ratio Z gauss (A + l)Z gauss (A — l)/Z gauss (A) 2 are linear functions of the matrix size A. This is the first 

^ee, however, |g] where the weak-coupling expansion in gauge theory is compared to the non-perturbative expansion in string 
theory over the world-sheet instantons. 
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sign of appearance in the context of matrix models of the Toda lattice equations, to be considered below. The 
quasiclassical expansion 1(1.9(1 of (|1.10() reads 

F gaues = ^ (| (logic - + £to*?) + 0(h°) (1.11) 

where the first logarithmic term in the r.h.s. comes from carefully treated large N asymptotic of the volume 
factor Vm- Expression in brackets in the r.h.s. of (|1.11|) . proportional to 1/ti 2 , is a "mostly trivialized" example 
of quasiclassical tau-functions to be in the center of our interest below. 



2 Matrix models and integrability 

Before discussing the quasiclassical limit of the matrix integral ((1.3|) . let us remind some its general properties. 
The A r2 -tuple integral can be computed, first using gauge symmetry, being reduced to the only TV-tuple eigen- 
value integral. The key properties of the last one are in fact encoded into a single (contour) integral, due to the 
fact that such integrals are directly related with the integrable systems of Toda family. 

2.1 Matrix ensembles and Toda lattice 

Let us start here will the well-known and widely used statement of [5], that matrix model partition (or gener- 
ating) functions satisfy equations of the Toda hierarchy, or, more strictly, that certain integrals over random 
matrices are particular tau-functions of the Toda lattice hierarchy. Using the gauge invariance $ — > of 
the "single-trace" potential 1)1.4(1 of general form, one can rewrite formula ((1.3(1 as an eigenvalue integral (in 
this section we put h= 1), i.e. the partition function ((1.3|) equals Z = rjy(i), where 

i r N 

^) = M E{d^- Wi ^)^W (2-1) 

l — l 

and 

A(0)=n(^-^) = det||^r 1 |l (2.2) 

i<j 

is the Van-der-Monde determinant. 

In order to verify ((2.1|) one has to use the gauge transformation, which allows to diagonalize the matrix 
$ by $ = W(j>U, where U £ U(N)/U(1) N belongs to some "orbit" of the unitary group U(N) and 4> = 
diag(</>i, . . . , 0jv) is a diagonal matrix. The Van-der-Monde determinant comes into the integration measure of 
((2.1(1 as a Jacobian of transformation 

**«l*=0 = {u*<t>u - 4>) l0 u= = n) -*P fi > 4>]ij = -itfUjfa - (f>j) (2.3) 

Analogous to ((2.1(1 eigenvalue representations exist also in the "multimatrix" case, where the integral is taken 
over some ensemble of random matrices with special single-trace potentials. For example, in the two-matrix 
model (the "complex-conjugated" version to be mostly considered below) one introduces two matrices $ and 
& with the interaction 

V($,&) = &$-W{$;t)-W(&-,t) (2.4) 

and the matrix integrals reduce to the integrals over their eigenvalues J d$d& exp (— V($>, $^)) = Tjv(i, t), 
where similar to ((2.1(1 

. N 

rw(t,i) = ^ / ]J (dV- y <«'*>) |A(z)| 2 (2.5) 

7=1 
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due to the Harish-Chandra-Itzykson-Zuber formula. In the representations and (|2.5[l the JV 2 -tuple matrix 
integrals are reduced to the TV-tuple eigenvalue integrals (in the last case the integration over the eigenvalues 
(zi, Zi) S C C C 2 is taken generally over some submanifold of C 2 of real dimension two), and this is a particular 
case of quite common in gauge theories localization procedure. 

The most effective way to calculate the eigenvalue integral is to use the method of orthogonal polynomials. 
The Van-der-Monde determinant (|2.2(l can be equally written as 

A(0) = detll^r 1 !! =det||P i _ 1 (&)|| (2.6) 
with any polynomials, normalized with unit higher coefficient 

P k (<P) = <p k + 0{<j> k - 1 ) (2.7) 

and for the computation of matrix integral l|2.I(l it is convenient to choose particular Pj (</)) being orthogonal 
polynomials with the weight, determined by the matrix model potential W(<i>) 

= ( Af>e- w ^Pi{4>)PM = %e 9l(t) (2-8) 



Substitute the Van-der-Monde determinant in the form l|2.6|l . where the polynomials satisfy (|2.8|l . into i|2.1[l . 
After expanding the determinant and using orthogonality l|2.8l) . the partition function l|1.3|l acquires the form 
of 

N-l 

r N (t) = TT e«<*> =detJTy(t) (2.9) 
to 

a product of their norms (|2.8|l or the determinant of the so called "moment matrix" 

H l3 {t)= f 4> l+ Je- w{ -^dcl> (2.10) 



Formulas (|2.9|) . I|2.10|l claim that matrix model partition function (|2.1|l . since it can be presented in the 

form of determinant of an operator, exponentially depending on times or parameters of the potential (|1.4|) . is 
a tau-function of Toda hierarchy and satisfies the corresponding nonlinear differential equations w.r.t. these 
times. 

To derive these equations, consider the Lax operator for the Toda chain, which in the basis of the orthogonal 
polynomials (|2.8|l is presented in the common form of a tri-linear matrix [5] 

0P,(0) = Pi+i(0) - Pi(t)Pi(cj>) + RityPi-ttf) (2.11) 

The first from the set of Lax equations dt k L — [L, 1Z o L fc ] , with k = 1, immediately gives rise to 

d 2 q t 



(2.12) 



the equation of motion of exponentially interacting particles, (with co-ordinates qt and momenta pi), or the first 
equation in the Toda-chain hierarchy 

Indeed, from (|2.8() . I|2.11|l one can easily establish the relations among the functions qi{t) and the matrix 
elements Ri(t) and Pi(t) of the Lax operator. First, 

- 1) = e 9l(t) = Ritye*- 1 ® (2.13) 

gives rise to 

R t (t) = e <r*(*)-9«-i(t) Q.IA) 
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Differentiating 1|2.8|) for i = j one gets 

± m _ «.§* . / ^ «" + 2P,^) . p,e» (2.15) 

where the second term in brackets disappears due to orthogonality condition and condition (|2.7(l . Now, from 
(|2. 15|) it follows that 

ft(t) = (2-16) 

is just the momentum of i-th Toda particle with the coordinate qi(t). Differentiating now 1|2.8[) with i > j, we 
obtain 

= j #e-£ ^* + P^\ (2.17) 

and using l|2.8|) for comparing this for j < i — 1 with (|2.11|) one gets 

BP 

5? = fl^i-i (2-18) 

OTi 



Now we are ready to differentiate (I2.11|) : 

, ap op +1 d Pi 



Pi + O^- 1 ) (2.19) 



dti dt\ dt\ 

Multiplying l|2.19[l by Pi integrating and using l|2.8|) and l|2.18|l one finally obtains the equation 

^ = R l+l - Ri (2.20) 
at i 

or, using Ij2.16|l and 1|2.14(1 . arrives to the Toda chain equation of motion (|2.12(1 . 

For the partition function of matrix model (|1.3fl . (|2.1() or tau- function of the Toda chain (|2.9() . equation of 
motion (|2.12|l can be rewritten in the form of one from the infinite set of the Hirota bilinear equations 

atf logrjvW = — — ( } 

In the two-matrix case instead of 12.8fl one should introduce the bi-orthogonal polynomials 

d 2 ze- v ^^P t {z)P,{z) = <%e« ( *' t_) (2.22) 

giving rise in a similar way (see [H]) to the two-dimensional Toda lattice hierarchy with the set of times, 
corresponding to the potential l|2.4(l . In particular, the second time-derivative Jp- is substituted by the two- 
dimensional Laplacian , and the first Hirota equation (|2.21|) is, correspondingly, replaced by 



a , r log T N {t,t) = — (2.23) 

2.2 Matrix integrals in Miwa variables 

An interesting phenomenon occurs, when we rewrite the matrix integral (|1.3J) in terms of the Miwa variables 

ED 

i fc = 4 0) -^X>7 fc -4 0) ^TrM- (2.24) 
i=i 
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or parameterize the times by eigenvalues of auxiliary L x L matrix M. For the background values let us take 
tfc = 5<5fc,2, i.e. consider the matrix integral around the gaussian potential i|1.10[l . 

Upon substitution of (|2.24l) . the eigenvalue integral l|2.1[) acquires the form 



w (o) ) = ^ /n ( 



i( > 3 wni-? = 



i r N f 



i,I 

A(<M 



/./,/ 



(2.25) 



where A(/x) = Hj< ~ ^j) an d A(</>, /x) are the corresponding L x L and (AT + £) x (JV + L) Van-der-Monde 
determinants. As in previous section it is again convenient to introduce the orthogonal polynomials Q2.8[l . I|2.7|l . 
which allow to rewrite (|2.25|l as 



Z(M,t (Q) ) 



x det ||Pj-_i( 



N 

no 



ie - w ^* m A x 



det 

(i,i),(j,J) 



Pj-1+N(<f>i) ■■■ Pj-l+N^l) 



(2.26) 



Expanding the determinants oi N x N and (N + L) x (N + L) matrices in (|2.26|l and using, as before, the 
orthogonality relations one gets 



r N {t) = Z(M,i<°>) = ^T!} \ d r ct \\Pj- 1+N ( N )\\ 



or, introducing ipj(ji) = fi N Pj +N (fJt), 



r N (t) _ det/j HV'j-i^i 



A(M) 



(2.27) 



(2.28) 



which is a particular well-known representation for the Kadomtsev-Petviashvili (KP) tau- function in Miwa 
variables (|2.24l) . For the gaussian case t^' = |#fe,2 the orthogonal polynomials can be identified with the 
Hermite polynomials 



rl)j(jn) = fi~ N / dx e - x *+» x x N+J 



fl — >-oo 



(2.29) 



which allows to rewrite l|2.28|) in terms of the matrix integral in external matrix field of general structure 

Z(A) OC fdX e -^(X)+TrAX (2 30) 



where for particular case of the gaussian model (|2.25(1 , as follows from the integral representation Q2.29jl , A = M 
and V(X) = \X 2 — NlogX. Existence of dual matrix integral representation (|2.3U|) follows from another 
integrable structure |13| . arising for the matrix integrals in external matrix field. 



2.3 Kontsevich integrals, free fermions and topological strings 

Let us now consider in detail the matrix model (|2.3l)|l . where the integral is taken over the L x L matrices, and 
the partition function is at the moment written up to a simple [M -dependent) normalizing factor. Similarly to 
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the expansion around the gaussian point in the previous section, the matrix integral 1)2.31))) can be considered 
as a (formal) expansion around its "classical trajectory" X = M , determined by the equation of motion 

A = V'(M) (2.31) 

The integral (|2~3n)l can be then expanded in the inverse powers of matrix M, i.e. the (gauge-invariant) partition 
function 1)2.30)) becomes the function of the variables 1)2.241) . and normalization is chosen in such a way, that 
this expansion starts from unity. 

In general case the matrix integral l|2.3U|l can be written as well in the form of determinant formula, like 
(|2.28|l . which states that it equals to a tau- function of KP hierarchy (more strictly for a polynomial potential 
V(X) of degree p it is the p-th Gclfand-Dickey reduction of the KP hierarchy). The determinant formula follows 
directly ^2] from evaluation of the integral 12.30fl 



AY p-Tr(VPO-TrAX) - 

A(A) 



rl W A (|)ll/ dx ie - v ^ +x ^ = A- 1 (A)dct||vI/ / _ 1 (A J )|| 



(2.32) 



where the first equality is a consequence of the Harish-Chandra-Itzykson-Zuber formula, and the elements of 
the matrix in the r.h.s. are expressed through the contour integrals (cf. with its particular case 1)2.2911 ) 

*,(A) = J dx z'e- v W+ A * = (J^j 1 * (A). (2-33) 

Putting on classical solution (|2.31|) A = V'(M), denoting the eigenvalues of M as {hi} and restoring the 
normalization, one gets the final result 

7(M\ « Tjyft) _ det fJ $/-i(^j) 

z{M) tiM^W)~ W) (2 - 34) 

with 

$ J (/i) = (V"(//)) 1/2 e v W-^'W* 7 (V'(//)) f/(l + o(-)) (2.35) 

which satisfies the Hirota equations of KP hierarchy, due to simple bilinear relations for the determinants, see 
details in |13j . The determinant formula 1)2.34)1 holds for any matrix size L and in this sense, unlike the matrix 
integrals ()1.3)) . H2.ll) . I)2.5|) . does not "feel" the L — ► oo limit of large size of the matrices we exploit intensively 
below. We immediately see, that representation ()2.28|) is a particular case of the formula 1)2. 34)1 . provided ()2.33)) 
are given by particular expressions ()2.29)l . 

The formulas of this section require few necessary comments (see |13) for their detailed discussion): 

• The (dispersionful!) Hirota bilinear equations, satisfied by the determinant formula ()2.34)) are in fact 
equivalent to the free field representation of the matrix elements 1)2.35)1 

$/(/!) cx (J|V>(/i)&|0) (2.36) 

with any operator of the form $ = exp 3 ] Aijipiipj^j . The determinant formula 1)2.34)1 then simply 
follows from the Wick theorem. 
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The particular functions (|2.35|) . possessing simple integral representation l|2.33[) are specified by string 
equation, which can be written in the form Z = 0, for the first-order differential operator w.r.to the 
times of the hierarchy 

m_x-r r J i \JL + l^f 1 v"to)-v"(i*j)\_ 9 



C - 1 ^ 1 Tr {v"(M)M^) dt n + 2^ {v"(m)V"(»j) ' in-iij' ) dti ( 2 ' 3? ) 

acquiring the more common form of a particular Virasoro operator = C_J for monomial potential 
V(X) = 



x p 



p+i 



Generally the Virasoro constraints Z = for n > — f can be used as equations defining (together with 
some extra conditions) the partition function of certain topological string theory |15| . or (p, g)-minimal 
string theory. Equivalently, one may impose only the constraint (|2.37l) together with the full set of bilinear 
Hirota equations of (p-reduced) KP hierarchy . The planar or quasiclassical contribution to the string 
partition function is then described by dKP (or dispersionless KP) hierarchy, to be discussed in detail 
below. 

A particular (p, q) string model is then defined as a q-th critical point of the p-reduced KP hierarchy, which 
is characterized as solution to all above equations in terms of the expansion near the point tk = r^^fc.p+g 
in the space of times or couplings. 

Matrix integrals l|2.3U[l are often interpreted as partition functions of the topological Landau-Ginzburg 
theories of a single field with superpotential V'(X), interacting with two-dimensional topological gravity. 
The target-space for these theories is a sphere with a marked point, or just a complex plane, and it can be 
conveniently identified with the rational curve y = V'(x), defined by an analytic equation (|2.3I|1 on two 
complex variables (x,y) = (A, A). This interpretation was proposed among the first examples of defining 
the nonperturbative string theory in terms of a (quasiclassical) tau-function, associated with a complex 
manifold endowed with certain additional structure: two given functions, or, better, differentials dx and 
dy with fixed (here vanishing) periods. The potential of matrix model in external field l|2.30[l V = J" ydx 
is then an integral of generating differential for the Landau-Ginzburg theory. We will see below, that this 
geometric language is universal and effective description for the quasiclassics of all matrix integrals. 

Integral representation (|2.33|) can be then thought as a particular case of the duality transformation 

oc [ d f j,{x)^i(x)e^ ydx (2.38) 



between the {q,p) and (p, q) models. The dual wave functions are naturally expressed as functions of dual 
local co-ordinates, x and y correspondingly. On classical level this duality therefore looks like change of the 
local co-ordinate, and again, we will see this below as general phenomenon in quasiclassical formulation 
of the matrix models. 



3 Planar limit of the one-matrix model 

Let us now come back to the simplest matrix integral (|I.3|I and turn directly to the discussion of the planar 
limit N — ► oo, h — > with the new parameter, introduced in (|f .8(1 . Nh = to, being fixed. In the limit N — > oo 
from the whole topological expansion (|f .71) only the first term Fq = T with g = survives, corresponding to 
the summing over the planar fat graphs, which can be drawn on sphere or plane without self-intersections. 
Simultaneously in this limit, since h — > 0, the matrix integral l|I.3|l can be calculated quasiclassically, i.e. by 
the stationary phase method. The study of extrema directly in the integral over matrices l|I.3|) is senseless, due 
to the huge gauge degeneracy of the extremal configurations in the space of all matrices, and therefore we turn 
below to the study of quasiclassical properties of the eigenvalue integrals (|2.I() and (|2.5(l . 
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Figure 1: Cuts in the eigenvalue plane, and cycles on the re-sheets of the hyperelliptic Riemann surface (|3 . 1 U|) of 
the one-matrix model. The A-cycles are conventionally chosen to surround the cuts or eigenvalue supports (on 
this picture supposed to be a segments of real line Imx = 0), while the B-cycles connect two different eigenvalue 
supports. 



3.1 The eigenvalue integral and complex curve 

The quasiclassical value of the eigenvalue integral 12.1fl at h — ► 0, can be found by studying the extrema of the 
effective potential 

W eB ((!>)=W(4>)-2hlogA(<t>) (3.1) 
or, more strictly, solving the stationarity equation 



(3.2) 

If ft = the interaction in the r.h.s. of (|3.2(l is switched off, and all eigenvalues 4>j are somehow distributed 
over the minima, or, in the most general setup, all extrema where W'(<f>) = 0. More convenient at TV — > co is 
to introduce the eigenvalue density 

N 

p{x) = hY J S{x-(t> j ) (3.3) 
j=i 

or the resolvent 

GW^Av-LA .[Em (3.4, 
\ x-&/ 9 Jc x-<f> 

defined on the x-plane with cutoff eigenvalue support, or, as we see below on the double cover of some complex 
a;- plane, cut along some segments C = (J . Cj, (see fig. 0), and normalized as 

/ dxG(x) = [ d<f>p(<j>) = t Q (3.5) 
2?" Jc Jc 

In these terms equation (|3.2|) can be rewritten as integral equation 

W'(x) = 2 / p ^ d f = G(x+) + G(g_), xeVacC (3.6) 
J c x - <f> 

where x± = x ± zO are two "close" points on two different sides of the cut - above and below if distribution of 
eigenvalues is real and cut is stretched along the real axis. This equation holds at any point of the eigenvalue 
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support C, consisting of several disjoint pieces for generic polynomial potentials. The formula l|3.ti[) can be 
further rewritten as an algebraic equation on the resolvent G, see e.g. ^| El El 120] 

G 2 - W'(x)G = f(x) (3.7) 

where fix) in the r.h.s. is a polynomial of the power one less than that of W'(x). Indeed, it follows from (|3.6|) 
that 

(G 2 - W'G)\ x=x+ - (G 2 - W'G) \ x=x _ = (G(x + ) - G(x-)) (G(x+) + G{x.) - W'(x)) = (3.8) 

i.e. the l.h.s. of (|3.7|) is a single- valued function f(x), whose properties (a polynomial of certain power) are 
determined from behavior at x — * oo. 

Equation l|3.7|l defines a hyperelliptic curve and, as we see below, the quasiclassical free energy T — F (or 
the first term in l|1.2|) 'l can be entirely defined in terms of the curve (|3.7|l and the generating differential Gdx, 
determined by the resolvent l|3.4|l and measuring the eigenvalue distribution. For the polynomial potentials 
equation (|3.7(l describes an algebraic curve (of finite genus n for the potential W of degree n + 1) and defines 
resolvent G as an algebraic function on the curve, which is double-cover of x-plane 2 . 

Introducing new variable by 

y = W'(x) - 2G (3.9) 

one can rewrite (|3.7|l as 

y 2 = W'{xf +4f(x) = R(x) (3.10) 

where R(x) is a polynomial of the degree 2n for the potential W(x) of degree n + 1. For such polynomial 
equation (|3.10(l is a canonical representation for the hyperelliptic curve of genus g = n — 1. 



3.2 Free energy: geometric definition 



By computation of free energy for the matrix model l|1.3(l in large N or quasiclassical limit, one usually under- 
stands solution to the variational problem 1)3. 2f> for the functional 3 



T = 



J dxp(x)W(x) - J dxi J dx 2 p{x 1 ) \og\x! - x 2 \ p{x 2 ) +^2u a (^J dxp(x)^j - S a 



(3-11) 



In addition to formal definition 1)2.1(1 we have assumed here some fixed distribution of all N — + oo eigenvalues into 
certain (finite number of) groups with the fractions ^2 a N a — N, implying an extra linear condition imposed 
on new quasiclassical variables S a = hN = to. Let us discuss, first, the nature and meaning of these groups 
and corresponding fractions. 

Intuitively it is clear, that in quasiclassical limit the eigenvalues form some condensates, which can be located 
around any minimum (or even near any extremum W'((f>) = as we determine immediately). If one turns off 
the Coulomb repulsion, coming from the Van-der-Monde determinant, by h — > the eigenvalues will be located 
exactly at these extrema. Therefore, for some small eigenvalue fractions this picture is evident, for the large 



oo it becomes a curve of infinite genus and nothing can be 



2 However, if one allows all possible long operators Tr^?^ for L 
said about the resolvent G immediately. 

3 By writing modulus in the argument of logarithms we generally imply choosing appropriate branch of the logarithm function, 
it literally coincides with modulus when all eigenvalues are supposed to be real. 
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A Im (j) 



Re (f 



Figure 2: Complex geometry of the potential W{(f>) = ^ t(f>. The level lines, passing through the critical 

points (f)± = ±\/i are Im</> = and (Re0) 2 — i = 0. The Coulomb repulsion stretches the eigenvalues 

along the segments of the level lines near the critical points. 



eigenvalue fractions it can be more complicated (the condensates can grow collide) and we will discuss some 
details of this below. The jumps between different condensates are exponentially suppressed, since they are of 
the order e - const / h ^ e -const-A^ anc j therefore, it is natural to introduce extra (to the coefficients of the potential 
1)1. 4p ) variables, or to fix the fractions of eigenvalues 

S a = hN a = [ dxp(x) = — <£ Gdx (3.12) 

JC a 27TI J Ca 

at each extremum, when solving the variational problem for the functional 1|3.11[) . say, introducing the corre- 
sponding Lagrange multipliers II a . 

Moreover, it is clear, that at switched off Coulomb interaction, the stationary configurations, in general 
position, consist of separated points - the extrema W'((j)) — of the potential, for large ./V each taken with 
certain multiplicity. However, the interaction smoothes the picture - instead of separate points one now gets 
the set of one-dimensional lines, to be possibly thought of as analogs of mechanical trajectories or Ising strings. 
For analytic continuation of the partition function in the space of parameters and deformation of the integration 
contours, this leads to appearance of the complex curves of statistical matrix ensembles. 

To understand this better, one should remember that as any integral, the eigenvalue integral l|2.1|l depends 
on the choice of the integration contours. For the potential of degree n + 1 the basis in linear space of all 
possible contours, when integral over each eigenvalue converges, is n-dimensional. It is easy to see, that each 
"steepest descent" contour can be drawn in the complex plane exactly passing through one of the n extrema 
points W'((/)) = 0, the recent detailed discussion of this issue can be found in 25,. 

For the simplest non gaussian potential W{4>) = ^ — t(f> this geometry is illustrated at fig. [21 Two critical 
points of potential are at <p± = iyt (at the figure we have chosen t to be real and positive). The steepest descent 
lines of ReVF, passing through the critical points </>±, are the level lines of ImW — Im</> ((Re<^) 2 — P m $) t) = 
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0. The integration contours, for which the eigenvalue integrals converge should start and end at one of the n+1 
sectors of the complex ^-plane, where KeW > (n = 2 for fig. [21 i.e. for this case there are three such sectors). 
Since integrands are holomorphic and the value of the integral does not depend on the local deformations of 
the contour, the steepest descent method "chooses" the integration contour from one allowed sector to another, 
passing through a critical point (two such possible contours are shown by dashed lines at fig- ■ As already 
mentioned above, the contribution of the Van-der-Monde determinant to the effective potential (|3.1[l leads to 
repulsion of the eigenvalues, so that from their locations exactly at the extrema W'(x) = they stretch along 
the segments of the allowed contours, see fig. 

The quasiclassical free energy (|3.11|1 is calculated at critical densities, i.e. when 



5T 
8p{x) 



W{x) + 2 / d4>p{4>) log \x - <j>\ - n Q = 0, xeC a (3.13) 
ic 



The right equality can used for computation of the Lagrange multipliers IIo,, which will be intensively used 
below, since they are equal to the partial derivatives of the free energy i|3.11fl w.r.t. the eigenvalue fractions 
j3~T2fl 

-IL a + dx——-^ r -!-=n a (3.14) 



dS a J c 5p{x) dS a 

since the variational derivative vanishes due to (|3.13() , which simultaneously determines the Lagrange multiplier 
itself 4 . 

Formulas i|3.11[) can be now accepted as a separate geometric definition of the planar free energy T . To make 
it strict, the solution to l|3.13[l . I|3.11[l should be supplemented by boundary condition - a choice of support C 
for the nonvanishing density p(x) ^ or non-analyticity of the resolvent G(x)dx (13. 4|) . This support (the set of 
" cuts" for the one- matrix model or " drops" in two-dimensional case of (|2.5|l , see fig. ^ and fig. [2J is introduced 
"by hands", but its length and form (or size and shape in two-dimensional case, to be discussed in detail below) 
is then determined dynamically. The density of eigenvalues in one-dimensional domain (union of "cuts") is 
expressed by (|3.7|l through the parameters of the potential and coefficients of auxiliary function 

n-l 

/(*) = E f* xk ( 3 - 15 ) 

fe=0 

related to the filling fractions. Note immediately that ambiguity in the choice of extra free parameters of the 
polynomial (|3.15|l . arising in the r.h.s. of l|3.7|l is exactly "eaten" introducing extra parameters - the filling 
fractions (|3.12() . 

Below we are going to show that the geometric definition of free energy (|3.14|l coincides with a particular 
case of so called prepotential or quasiclasical tau-function, which can be defined in principle for generic complex 
manifold endowed with some extra structure, usually encoded into a meromorphic differential form. For the 
planar limit of one matrix-model discussed in this section the complex manifold can be identified with the 
hyperelliptic curve (|3.7|) . H3.10jl . 



Example: the gaussian case 

Let us now solve explicitly the above equation ({3.6(1 for the simplest gaussian potential W — \x 2 . The r.h.s. 
of l|3.7|) or the polynomial H3.15fl in this case (ra = 1) is some constant, and the solution to 13.71) reads 

G=| + i v / ^+47^| + ^^ J ^ (3 . 16) 

4 Note also, that the Lagrange multipliers in 13.111 are also necessary for satisfying dispersionless Hirota equations and their 
analogs, i.e. they play the same role as careful normalization in the measure 11.51 dictated by 12.211 . H2.23I . 
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i.e. it has a single cut in the eigenvalue plane, with the only parameter a 2 = —4/, identified with the squared 
half-length of a symmetric around x = cut in the x-plane. It can be determined from the normalization l|3.5|l . 
which gives rise to a 2 = 4io- In this case the matrix model curve 1(3. 7|) . 1(3. 1U|) is rational or has genus g = 0, 
and the density equals 

p(x) = ^-^a 2 -x 2 (3.17) 
and one can check that the partition function for this case is expressed by the following formula 



tl A . 3 



ii=0 



= I (logic- 5) CUS, 



coinciding with the first term of the quasiclassical expansion of (|1.10|) at t\ — 0, corresponding to symmetrically 
located cut C = [—a, a] in the £-plane and the Wigner eigenvalue density 1(3.17(1 . Formula 1(3.18(1 follows from 
the only in the gaussian case relation 1)3.14(1 expressing derivative of free energy w.r.t. the total number of 
eigenvalues S = to 



-W(x*) + 2 / d(f>p((f>) \og\x* - (j>\ = -/ \Ja 2 - x 2 log xdx = t (logt - 1) (3.19) 

OTO Jc JO 4a 2 =t 

Integration of l(3.19|l leads to 1(3.18(1 . where the dependence on t\ is trivially restored by the shift argument (e.g. 

___ . . . _*0*1 

from 1(1.10(1 it immediately follows, that -Z gall ss = e ^ ^ g auss| tl=0 - 



4 Quasiclassical tau-function 

The quasiclassical free energy of the matrix model, defined by formulas 1(3.14(1 is a particular case of so called 
prepotentials, or more strictly, quasiclassical tau-functions, introduced in |2.'ij . In order to continue discussion 
of the geometric properties of matrix models we need to learn more about the properties of quasiclassical tau 
functions and to remind, first, some necessary facts from geometry of complex curves. 



4.1 Riemann surfaces, holomorphic and meromorphic differentials 

Let us start with reviewing some well-known facts from the theory of Riemann surfaces or one-dimensional 
complex manifolds - the complex curves. The topology of compact oriented Riemann surface is characterized by 
a single non-negative integer - genus g, generally seen as number of "handles" , attached to sphere. The simplest 
examples are sphere itself with no handles or g = and torus with g — 1 , an example of Riemann surface with 
g = 3 can be found at fig. [21 

On genus g Riemann surface S s there exists 2g independent noncontractable contours which can be split 
into pairs, usually called by A = {A a } and B = {B a } cycles, (see fig. El a = 1, . . . ,g, with the intersection 
form A a o Bp — S a p. The basis in the space of of one- forms is naturally dual to the basis of one-dimensional 
cycles, e.g. the holomorphic (or the first kind Abelian) differentials d(duj a ) = are canonically chosen to be 
normalized to the ^-cycles 

f duj a = S a p (4.1) 

JAp 

Their integrals along the B-cycles form then the period matrix 

duj a = T a p (4.2) 

B/3 
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Figure 3: Compact Riemann surface of genus g — 3 with chosen basis of A and B cycles, with the intersection 
form A a o Bp = 8 a p. Relabelling of A- and £?-cycles is called a duality transformation. 




Figure 4: Cut Riemann surface with the boundary <9E. The integral over the boundary can be divided into 
several contributions (see formula l|4.34|l ). In the process of computation we use the fact that the boundary 
values of Abelian integrals on two boundaries of the cut differ from each other by the period integral of the 
corresponding differential dio a over the dual cycle. 



(a symmetric g x g matrix with the positive definite imaginary part), which is one of the most important 
characteristics of complex curve. The symmetricity of period matrix follows from the Riemann bilinear relations 
for holomorphic differentials 

0=1 dup A du) 1 = I (j> duip 6 duj-f — (j> dto 1 (j> dup J = T ( g 7 — T^p (4-3) 



s 

coming from the Stokes theorem on the cut Riemann surface, (see fig. 0}, and vanishing of any (2,0)-form, in 
particular dtop A dto 7 , on a one-dimensional complex manifold. Analogously, the imaginary part of the period 
matrix is positive-definite, since lmT a p cx J s du> a A dwp. 

In contrast to holomorphic, the meromorphic differentials dQ are analytic everywhere, except for the finite 
number of points P, where they can have poles and BdQ, cx S(P). Their canonical form is fixed not only by 
half of the periods (as for the holomorphic ones, see l|4.1(l '). but as well by their behavior ("main part") at 
singular points. The common classification for the meromorphic differentials includes the second kind Abelian 
differentials, with a single singularity at some point Pq of the form dilk ~ ^+1 + • • •> k > 1, where by dots 
we denoted a nonsingular part, and ( is a local co-ordinate near the point Po : C(-Po) = 0. Note, that this set 
starts from differentials with the second order poles, since there is no meromorphic differential with a single 
first-order pole, due to vanishing of the total residue. Instead, for taking into account the first-order poles, it is 
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convenient to introduce the third-kind Abelian differentials dQ±, with two first order poles at some points P± 
and, opposite by sign, unit residues there. 

Note now, that to each so defined second and/or third-kind Abelian differential one can add any linear 
combination of the holomorphic ones, without changing their behavior at the singular points. To fix this 
ambiguity we should impose g constraints (the dimension of the space of holomorphic differentials) to their 
periods, the canonical choice of such constraint is vanishing of all their A-periods dfl = 0. This means, that 
as a strict definitions of the Abelian differentials one should write 

d^k p + . . . , / dfl k = 0, k > 1 (4.4) 



A 



for the second kind, and 



dn ~ ±^ + ..., rfSlo = 0, (4.5) 



A 



for the third kind, where C±(P±) = 0. The bipole differential (|4.5|) can be also presented as 

dna=diog wM (4 - 6) 

where E(P,P') is the so called prime form on £ x E, the — i-bidifferential with the only zero at coinciding 
arguments, see details and its properties in the book [22] . 

In the same way as with (|4.3|) . one can write down the Riemann bilinear relations for generally meromorphic 
Abelian differentials. For example, for the first and third kind Abelian differentials one gets 







/ (Lop A dflo = ( f dujp d> dflo — (p dQ ® dui^ 
Js Z \J A a JB a J A a JB a 



+res oo (duj0) j dtto — Tes 00 (dflo) / diop = d> dflo — j dup 



(4.7) 



which is widely used the theory of Toda integrable systems, appearing rather naturally in the context of matrix 
models. 

In practice, complex manifolds can be effectively described by systems of polynomial equations in multi- 
dimensional complex Euclidean spaces C K , to get a complex curve one should have K — 1 equations for K 
variables. The coefficients of these equations play the role of moduli of the complex structures on a given 
Riemann surface, for a genus g > 1 surface the complex dimension of moduli space is dime M. g = 3g — 3, that 
can be found from the Ricmann-Roch theorem. The Riemann surfaces with different complex structures cannot 
be holomorphically mapped to each other, i.e. are different complex manifolds, though can be isomorphic as 
real surfaces. Note also, that it means that matrix elements of the period matrices (|4.2|) are not functionally 
independent, since there are totally g{g + l)/2 functions of only 3g — 3 moduli parameters, except for the cases 
g = 2, 3 being solutions to the quadratic equation g(g + l)/2 = Zg — 3 and exceptional case 3=1 with the only 
modulus of a complex torus. 



4.2 Definition of quasiclassical tau-function 

Let S be a Riemann surface of some finite genus g 7 and suppose we have defined on S two meromorphic 
differentials, say denoted dx and dy, with the fixed periods. As we already discussed above, without loosing 
of generality one can always take the canonically normalized meromorphic differentials, with, say, vanishing 
A-periods. However, if we simultaneously fix the rest -B-periods, this gives already nontrivial constraints on the 
moduli of S. In practice this restricts us onto some subspace in the total moduli space, and the dimension of this 
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subspace can be estimated in the following way: the differentials dx and dy are only denned up to multiples, and 
fixing the B-periods of dx gives g — 1 constraints (taking account of the scaling freedom) , while the -B-periods 
of dy give a further g — 2 constraints. Altogether we come to a system with (3g — 3) — (g — 1) — (g — 2) = g 
parameters, i.e. the number of moduli equal to the genus of E. 

A complex curve endowed with two meromorphic differentials with the fixed periods, or, briefly, the g- 
parametric family of curves defines an integrable system 5 . In practice, one considers a g-parametric family 
of Riemann surfaces E, endowed with a generating meromorphic differential, given, up to a constant, by the 
formula 

dS ex ydx, SmoduiidS — holomorphic (4-8) 

where y(P) = J P dy, P G E is an Abelian integral on E, and by variation in moduli we mean any local variation 
of the rest g parameters of the family, which "survives" after one fixed all periods of dx and dy. 

A general definition for the prepotential of an integrable system for a (^-dimensional family of smooth curves 
E 9 of genus g endowed with meromorphic one- forms l|4.8fl . is of the following form 

S = I dS, 

dT ^ , (4-9) 



as 2m V s 

where A and B are dual cycles in the homology basis. The prepotential is defined by ()4.9|1 locally on the 
moduli space of E, or better, on the Teichmullcr space with the fixed basis in i? 1 (E,Z). For complex curves, 
where it is natural to add to 1)4. 9f) the variables associated with the degenerate cycles or marked points (as 
in the simplest example of gaussian matrix model l|3.18[) we have considered above), the exact definition of so 
generalized prepotential was given in , where it was called the (logarithm of) quasiclassical tau- function or 
tau-function of generalized Whitham hierarchy. The consistency of (|4.9[1 is guaranteed by integrability condition 

2mT afl (4.10) 



dSadSp 

following directly from (|4.8|l and (|4.3|) or symmetricity of the period matrix T a p of the Riemann surface E s . 

For degenerate curves or even the curves with the marked points, the formulas (|4.9|l should be treated more 
carefully and cannot be applied directly. A degenerate ^4-cycle turns into a pair of marked points P± (a solitonic 
limit of a handle) and for such degenerate handle formulas 1)4. 9JI turn into 

<o = - — resp dS = — - — rcsp_dS 
2m 2m 

dT_ 

dt a 



(4.11) 

Am I dS 



where the last expression is naively divergent and should be applied only after extra careful definition, where the 
divergent integral is replaced by appropriate finite quantity, like it happened in the simplest gaussian example 

In most general setup, one should complete the definition l|4.9|) by the time-variables associated with the 
second-kind Abelian differentials with singularities at a point Pq 

t k = — 5— res Po £~ k dS, k > 
dT (4-12) 

— = — res Po £ k dS, k > 
oth Zm 



5 In the most general spirit of the Liouville theorem, one can say that on g-dimensional subspace of moduli space of S one can 
always choose g independent functions - Hamiltonians or actions, while the co-ordinates on Jacobian Jac of E, a g-dimensional 
complex torus, play the role of complexified angle variables. 
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where £ is an inverse local co-ordinate at Pq: £(Pq) = 00 ■ The consistency condition for (|4.12l) is ensured by 



' hI \ves Po (edil n ) (4.13) 



dt n dtk 2ni 

and symmetricity of l|4.13|l is provided by fl n = ff, where operation + leaves only the main, singular at Po, 
part of the function In other words, the differential dQ is a meromorphic on E g with the only singularity at 
Po, i.e. is the second-kind Abelian differential J23J. 

Altogether, the basis of "flat" times (S,t,io) exactly corresponds to the basis of the first, second and third 
kind Abelian differentials. For the generating one- form Ij4.8|l one has 



dS a 
together with 



= duj ai a=l,...,g (4.14) 



<>llS dfl (4.15) 



dt 
and 

^ = dn kl k>\ (4.16) 

dt k 

which are "dual" formulas to the first lines of Q4.9p . (|4.11|) and (|4.12|) . and derivative over moduli is taken at 
constant local co-ordinate £ (that corresponds to choice of connection on moduli space of £). For the rational 
curve So and a single marked point Po the variables l|4.9|l and (|4.11|) are absent, and formulas (|4.12l) define the 
tau-function of dispersionlcss dKP hierarchy. In this case the singular at Po part £™ is a n-th degree polynomial 
in some global uniformizing co-ordinate A. 



4.3 Prepotential beyond one complex dimension 

The quasiclassical tau-function (|4.9|l can be considered as a particular case of so called prepotentials of complex 
manifolds, which exist as well for higher complex dimensions, like in the case of the Calabi-Yau three-folds, 
being one of the most interesting examples for the purposes of string theory. A general definition of prepotential 
contains: 

• A complex manifold E, homology basis with symplectic structure A a o Bp — <5 Q/ g, moduli space M. of 
complex structures on S; 

• Generating differential form fi; 

• The set of period variables S a — § A f2; 

• Dual variables - the dual periods II Q = § B f2; 

• Integr ability condition = §§^, following from bilinear relations J s 5Q A 5tt = 0; 

If so, like in the one-dimensional case considered in previous section, one can argue, that there exists a 
prepotential T: II Q = . As in the one-dimensional case, the whole picture is defined only locally on the Te- 
ichmiiller space, but on the moduli space of complex structures it is consistent with the duality transformations: 
A «-> B, S <-> II, in the sense of the Legendre transform of the prepotential 

F^F + ^Salla (4.17) 
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or, generally, the whole "electro-magnetic" duality of Sp(2h,Z), where 2h is the dimension of the homology 
group 

However, there are few very important distinctions to be necessarily commented. For the moduli spaces in 
the case of Calabi-Yau three-folds one has literally h = dimAl = /i^ 2,1 ) + 1 = #(A — cycles) coincidence of 
dimension of moduli space with the number of A-cycles, unlike the case of complex curves, where dim M = 3<?— 3, 
while #(A — cycles) = g and there is discrepancy with naive counting. Therefore, in the latter case one has 
necessarily to restrict the family of £ to some subspace of moduli space, the most invariant way to do this is, 
as already mentioned above, to consider a pair of meromorphic one- forms with the fixed periods, i.e., if on the 
curve £ 

3 f^i, 2 = {dX, dz} : <j> Vt\ t 2 = const 

this leaves exactly g moduli from 3g — 3 and this is the most "rough" definition of an integrable system - at 
the level of counting, see footnote [SJ However, for higher dimensional complex manifolds the number of cycles 
usually exceeds the number of moduli. 

For generating differential forms one gets the following picture 

• For the Calabi-Yau 3-folds ft S H 1 - 3 * * 1 is a unique (up to multiplication by a factor) holomorphic (3, 0)-form 
and relation J^, 5S1 A 5S1 = follows from the decomposition 5£l € H^ z ^ © H^ 2 ' 1 ^; 

• For the complex curves f2 = dS is a meromorphic (1, 0)-form, in addition it may have nontrivial residues 
res(dS'), res(£,~ k dS), which play the role of the "generalized" periods l|4.11[) . (|4.12l) . 

Hence, in the last case we restore the picture of the previous paragraph. 
4.4 Tau-function of the one-matrix model 

Now let us go back and continue discussion of the one-matrix model free energy, using the new notions introduced 
above 6 . The functions G and y (see Ij3.7l) and (|3.10(0 are defined on the double cover of the complex x-plane 
and on each sheet of this double cover they acquire discontinuities along the segments or cuts, proportional to 
the eigenvalue density p(x), vanishing outside the cuts (see fig. Pi . Therefore, the eigenvalue distribution (|3.12l) 
can be described by the periods of generating differential l|4.8[) 

S a = j dS (4.18) 

where the contour integrals are taken around the eigenvalue supports C a , to be identified (except for one of the 
supports, i.e. k = 1, . . . , n — 1) with the set of canonical A-cycles on the curve l|3.1(Jfl (see fig.^l. The generating 
differential Gdx with the resolvent l|3.4|l for the one-matrix model can be equally chosen as 

dS = — ydx (4.19) 

47T 

where now x and y are the co-ordinates on £ determined by l|3.10|l . since § Gdx = — \ § ydx, see l|3.9ll 7 . Both x 
and y are globally defined meromorphic functions on l|3.10|l , so that the integrals <f dx = and <f dy — obviously 
vanish for any choice of the closed contours. The relations (14. 14ft are certainly valid for the differential l|4.19|l . 
when the derivatives are taken at fixed residues and coefficients {ti} of the potential i|1.4|) . 

6 The geometric formulation for the free energy of one-matrix model was first proposed in 1181 and was discussed later in 
pll20ll2T) . 

7 The question about residues 14.111 . 14.121 is only slighty more delicate. For the symmetric choice of generating differential 
14.191 one can compute them at any of two infinities x = 00, while in the asymmetric case 13.41 they should be computed only at 
infinity on unphysical sheet, where the resolvent 13.41 is singular. 
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Figure 5: The calculation of the integral in (|4.22Jl . The integral over the cuts Cfc for k = 1,2 (fat lines), 
surrounded by the A/j-cycles, are replaced, first, by the integral over the B'-cycle around the logarithmic cut 
(the dashed line). The latter is transformed to the integral along the contour between the points x \ and x 2 
(without logarithm in the integrand), which is half of the period taken along the B-cycle, going one way between 
the cuts Ci and Q 2 on one sheet and vice versa on the other. 



Now, an important fact is that the Lagrangian multipliers (see (|3.13f) 'l can be also rewritten as period 
integrals of the same generating differential H4.19fl over the dual B-contours (see again fig. ^) 

n Q = -W(x* a ) + 2 / dxp(x) log (x - <) (4.20) 
Jc 

where the points x* a € C Q can be conveniently chosen, each on a-th piece of the support; then the modulus can 
be "forgotten" (see footnote |3J) and II Q become holomorphic quantities. 

Formula 14.20|) can be transformed easily to the contour integral (for simplicity, we consider explicitly the 
case of only two cuts, as on fig. the generalization is straightforward). In this case we have two fractions 
Si + S2 = to and it is more convenient to choose the independent variables to and S = S 2 . Then, ^2 a Tl a S a — 
(II2 — IIi) S + Ilito and variation of the free energy 1 If) w.r.t. S gives 



dT_ 

dS 



= W{xl)-W{xl)+2 / dxp(x) log -^4 (4.21) 
„._i ,Jc„ x — x l 



a=l,2 



Since the integral over the contour around total support vanishes j> c dxG{x) log x x _ x x l = (this contour can be 
deformed to the contour around the point x = 00 where the integrand has no residue), the last expression can 
be transformed to (see fig. EJ) 



/ dxp(x) log - — — = — !— <f> dxG(x) log - — — 
a ~i,2 Jc « x-x 2 ~^^ 2m Jc a x-x 2 

= L I dxG(x) log = dxG(x) 

2tti J b , x ~~ x 2 J x , 



(4.22) 



Altogether, this leads for (|4.21() to the formula 

W{x\) - W(x 2 ) - 2 / 1 dxG(x) = I ' 1 ydx = -\ I ydx = 2iri I dS (4.23) 
JxX Jxi Jb Jb 



dT_ 

dS 



where the contour B goes between x\ and x 2 and back on the other sheet of the curve (do not mix with the B'- 
contour around the points xl and x 2 \). For generic number of cuts one should similarly replace J2a=i ^-aS a = 
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(Ill — n„) Si + . . . + (II n _i — n„) S n _i + n„to and proceed to the trick with the logarithmic cuts. This leads 
to identifying free energy of the one-matrix model with the prepotetial l|4.9[) . 

To the set of parameters <|4.18|) one should also add 8 the total number of eigenvalues l|3.5fl 

to = ^re Sco (dS) = 2 f n - 1 (4.24) 

which is exactly the zero time (|4.11JI . and the parameters of the potential 1)1. 4JI . rewritten in the form of the 
times H4.12JI . i.e. as 

t k = ^ res °° (x~ k dS) , k = 1, . . . , n (4.25) 
since x plays the role of the inverse local co-ordinate at the infinity point x = oo. Then, explicitly 

71-2 



d n ddS 



dt 



x n x dx n sr^dfkX k dx , A „„. 

(n + l)t n+1 + 2 z2t^T (4 - 26) 



s,t v t^ dt o y 



fc=0 

and the dependence of {/&} with k = 0, 1, . . . , n — 2 on to is fixed by the period integrals 

x n ~ 1 dx n s—^dfkX k dx\ „ , , „„, 

(n + lt n+1 + 2^^ =0 4.27 

which for a = 1, . . . ,n — 1 gives exactly n — 1 relations on the derivatives of /o, /i, . . . , fn-2 w.r.t. to. 
For the derivatives w.r.t. parameters of the potential (|4.25|l . one gets 

W'{x)kx k ~ 1 dx ^? dfj x k dx 



dds 

d\lk = 



dt k 



2 i^w h ^r (4 - 28) 

s,t v j=0 ui k y 



analogously obeying 



dQk= r w^kx^dx^^r xHx^ Q 

JA a V dt k J Aa y 

and this is again a system of linear equations, resolved for To complete the setup one should also add to 

3T 1 

— = re Soo (x k dS) , k > (4.30) 

dt k 2ni 1 ; y J 

and the following formula 9 

dT f 

ss:= n "=L dS ' 4 ' 31 » 

with the appropriately chosen contour B' n passing from oo_ on the lower sheet to oo + on the upper sheet 
through the n-th cut (we again remind that, instead of £o, the parameter S n = to — Ylj—i S a can be used 
cquivalcntly). In the case of differential (14.191) for the one-matrix model these definitions requires some care 
due to divergences, however, in the geometric context of the AdS/CFT correspondence one considers an 
example of the generating differential with regular behavior in oo± , and the formulas like (|4.31l) can be literally 
applied. 

8 By the oo-point in what follows we call for short the point oo+ or x = oo on the "upper" sheet of hyperelliptic Riemann surface 
13.101 corresponding to the positive sign of the square root, i.e. to y = +\JW (x) 2 + 4f(x). 

9 Naively understood the integral in 14.311 is divergent and should be supplemented by a proper regularization. This subtlety 
can be just ignored, if one needs only the residue formulas for the third derivatives, (to be considered in the next section), say for 
the WDVV equations 1241 . The simplest way to avoid these complications is to think of the pair of marked points co and oo_ 
as of degenerate handle; then the residue 14.241 comes from degeneration of the extra A-period, while the integral 14.311 from 
degeneration of the extra B-period, see fig. Q For practical purposes this divergent period can be always replaced by a finite 
quantity, see e.g. the gaussian example 13.191 . 
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4.5 Residue formula 

Now let us discuss one of the most universal and nice formulas for the quasiclassical tau-function - the residue 
formula for its third derivatives 

d 3 T 1 f dHjdHjdH K \ 

dT ld TjdT K = 2^*=° { d X dy ) (4 ' 32) 

which is universal for dependence of T upon any variables {Tf }, with {dHi} being the corresponding differentials. 
Existence of simple formula (|4.32[1 for the third derivatives together with absence if similar expressions for the 
higher derivatives has definitely the string theory origin, coming from the world sheet theory, where the three- 
point correlation functions on sphere play a distinguished role. 



4.5.1 Holomorphic differentials 

Let us derive the formulas for the third derivatives of prepotential following the way proposed by Krichever 
in (231 7 an d presented explicitly in 24 . We first note, that the derivatives of the elements of the period 
matrix l|4.10|l can be expressed through the integral over the boundary <9E of the cut Riemann surface E (see 

fig. a 

-r-^- ee d 7 T a p = / dyduj a = - I ujpdjdua (4.33) 

U^l JBp JdT, 

where up — J F dwp are the Abelian integrals, whose values on two copies of cycles on the cut Riemann surface 
(see fig-EJ are denoted below as ujg. Indeed, the computation of the integral in the r.h.s. of 14. 3311 gives 



/ topd^dwa = ( / ujtd^duja - u)g d 1 dui a - / ujtd k duj a 

Jot, ~ \Jb p Jb p J - \Ja p 

>J i> (f> duup dyduja — 2_. f f dujp d, 
p Jb p \Ja p J p Ja p \Jb p J 



A„ 



dui a = (4.34) 

One can now rewrite formula (|4.33|) as 

d-yTap = — Ljpdjdu) a = / d-fUspdiVa = y_\ resdA=o (d^updwa) (4.35) 

J OS JdT, 

where the sum is taken over all residues of the integrand, i.e. over all residues of dyUip since the differentials du> a 
are holomorphic. In order to investigate these singularities and clarify the last equality in l|4.35|l . we discuss 
first the strick sense of the derivatives <9 7 w.r.t. moduli, or introduce the corresponding connection on moduli 
space of E. 

To this end, let us introduce a covariantly constant function of moduli, say the Abelian integral x — J P dx, 
i.e. choose such a connection on moduli space, that <9 7 x = 0. Roughly speaking, the role of covariantly constant 
function can be played by one of co-ordinates, in the simplest possible description of complex curve by a single 
equation on two complex variables - any of two functions constrained by the equation F(x, y) = 0. Then, using 
this equation, one may express the other co-ordinate y as a function of x and moduli (the coefficients of the 
equation). Any Abelian integral u>p can be then, in principle, expressed in terms of x, and in the vicinity of 
critical points {x a } where dx = (in general position) we get an expansion 

LU/3{x) = UJp a + Cp a \/X - X a + ■■■ (4.36) 
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whose derivatives 

d^p = d y uj p \ x=const = - < ^— d 1 x a + regular (4.37) 

possess the first order poles at x = x a , and they are written up to regular terms which do not contribute to the 
expression (I4.35|) . The exact coefficient in (|4.37H can be computed relating the function x with the generating 
differential dS = ydx. Using 

y(x) = T a y/x - x a + . . . (4.38) 



where T a = J]\ b + a {x a - x b ) or 



together with 



and 



^ V{X) = -2^TM +reSUkr (4 ' 39) 



p 

dy = a = dx + regular (4.40) 
2y x — x a 

dwp = Cl3a dx + ... (4.41) 
2^/x -x a y ' 



and, following from (|4.19|) and (14. 14ft expansion 

duj 1 = d 1 dS = - ^^ Xa = dx + regular (4.42) 
zy x x a 

one finally gets for i|4.35|) 



^ / c/3 a d 7 x a , \ sr^ f duj P 7 n A ( dio a duipduj. 

res (a^dw,.) = ^ res ( k« ) = L rcs {j^^x. j = ^ res [~^^- 



(4.43) 



In the case with hyperelliptic curves (|3. 10|) for one-matrix model, this general derivation of the residue formula 
can be replaces by based on the formula [2"2"j 



dT aP 
dx a 



L) a {x a )u)p{x a ) (AAA) 



where Cj a (x a ) 



duj a (x) 



d\/x — x a 



is the "value" of canonical differential at a critical point. 



4.5.2 Meromorphic differentials 

Almost in the same way the residue formula can be derived for the meromorphic differentials. One gets 

= ^re Soo (x k dS) , k > (4.45) 
ot k 2ni v ' 

therefore 

d 2 T I 1 

-tt—^— = 7r^ res °° (x k dn n ) = — res^ ((fl k ) + dQ n ) (4.46) 
otkOtn 2-ni v ' 2ni 

where (£lk) + is the singular part of the integrated one-form dQ k - Further 

^-resoo (x k dtt n ) = reSoo (x k j = -reSoc UdQ k ) + ^^- j = -reSoc (dVl k ^^- j (4.47) 
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The last expression can be rewritten as 

- res - fait) = L ( dQk ik) + ^ res - ( dnk ^t) = Eres ^ ( dQk ^t 

since f as = due to £^ dn n = 0, (cf. with 

I "4*- = ? (1 - 1 n »-5» - e (1 - L n '^ da - 

- E £ (1 «<) - E £„ (1 **) - 

= E (£ «.) £ - E *<) £ 4;*. ,^ = J=0 » 

Now, as in the holomorphic case, one takes the initial terms of the expansion 



An. = An- 
dt k 3 dt k 1 



x=const ~ Xa dt k 

Then, using ((OS)) . 1 (05) 1 and ((PU)) together with 



and the relation, following from Q4.16|) . 14.39|) 



one gets for l|4.48|l 

, w. a , = ^— Vres Ka d^fc^— = -— > rcs 2;a dil k - = — 

otkOt n ot m 2m ' \ af m / 27Ti ' \ 2ya; — cc a at 

1 x -> / d£l„ aa; a \ 1 x / d£lkd£l n dQ, rt 



27ri ^— J Xa \ dx dt m J 2m ^— ^ 21 a \ dxiiy 



(4.48) 



(4.49) 



n„(x) = Q na + jnaVx - x a + ... (4.50) 

and, therefore 



ry. Qx 

' a + regular (4.51) 



^ = vfe da;+ - (4 - 52) 



dil k = Ad5 = --A^= ^ + regular (4.53) 
aifc 2y/x - x a ot k 



(4.54) 



The derivation of the residue formula for the set of parameters including to corresponding to the third-kind 
Abelian differential i|4.26[) can be performed almost identically, and in a similar way, one proves the residue 
formula for the mixed derivatives. Thus, we finally conclude 

d 3 T _ 1 x ^ / dHidHjdHx 
dTjdTjdTK ~ 2^i^ Tehxa \ dx~chj 

Xa (4.55) 
1 v - ^ (<M>Ak a \ y^t,2, / s A , s,, \ \ - H I (x a )H J (x a )H K (x a ) 

= 2^^ [-^JWV = = E Ub ^ xa _ xb)2 

for the whole set of variables {!/} = {t k , to, Si} and corresponding to them one-forms {dHf\ = {dflk, dtto, duJi}. 
In the second line of l|4.55[) we have introduced the meromorphic functions 

*'<•>- < 4 - 56 ' 

for any (meromorphic or holomorphic) differential on hyperelliptic curve i|3.10[l dHj = Hi (x) ^ . 
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5 Complex curve of the two-matrix model 



Now let us turn to the case of two-matrix model (|2.5|) . As for the one- matrix model we first discuss the auxiliary 
complex manifold E, and then turn to the quasiclassical tau- function. 



5.1 Quasiclassics of the two-matrix model 

Consider the free energy of the two-matrix model (|2.4|) . (|2.5|l in the planar limit. Again, as in the one-matrix 
case in the planar or quasiclassical limit one can replace the direct computation of the eigenvalue integral (|2.5|> 
by solution of the corresponding variational problem. At least for the real problem in two-matrix model, with 
the mutually complex conjugated eigenvalues 10 the analysis becomes easier due to the basic property of the 
complex logarithm 

Alog|z-z'| = 2nS^ (z - z') (5.1) 

i.e. the two-dimensional Coulomb repulsion of the complex-conjugated eigenvalues is now governed by a kernel 
of the operator inverse to the two-dimensional Laplacian A = 4d z dz- Therefore for the potential 1)2.4(1 the 
eigenvalue density is some constant inside the support or (two-dimensional) eigenvalue domain and vanishing 
outside, i.e. 

, s f 1, z € D 
p(z,z)<xi u 5.2 

PV ' \ 0, z i D K ' 

and the corresponding real problem reduces to finding the potential energy of some domain (a set of "drops") 
in complex z-plane, filled by charged liquid with constant density of charge, see fig. [jjf . In contrast to one- 
dimensional case, now the eigenvalue density (|5.2|l is constant and the nontrivial information about their dis- 
tribution is encoded into the area and shape of domain, expressed through the coefficients of matrix model 
potential l|2.4|l (more strictly of its harmonic part) . 

The stationarity equation 

z l = W\z l )~hY J -— (5.3) 
following from the two-dimensional frcc-cncrgy functional, corresponding to the effective potential 

(5.4) 



1<J 

for the two-matrix case (|2.5(l can be re-written in the form 

z = W'{z) - G(z) (5.5) 

(together with the complex conjugated formula) after introducing, like in the one-matrix case I|3.4|l . the 

density l|5.2|l and the following function 

G(2) . fi ( Tr -LA . I / s<M)£i I f j£L m 



10 Sometimes it is also called, in a rather misleading way, the "normal" matrix model. We stress here, that studying the partition 
function 12.51 with the particular "cross-term" Tr (<J>t<j>J m t ne potential 12.41 the only essential thing is the choice of "real section" 
or the particular class of integration contours over each Zj and Zf integration variable. The difference between "normal" (with the 
commuting matrices $ and $') and usual two-matrix model arises only when studying multitrace correlators, which are beyond 
the scope of geometric description discussed here. 

llr This problem has attracted recently lots of attention in the context of so called Laplacian growth, see, for example 1271 and 
references therein. 
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Figure 6: Drops in the eigenvalue plane of the real problem in two-matrix model. In what follows this set of 
drops will be referred to as a domain D. 

or the (holomorphic) resolvent for the two-matrix model. When interaction in r.h.s. of 1)5. 3[) is switched off at 
h — ► 0, the resolvent vanishes, and l|5.5l) can be symmetrically written as 

(z-W'{z)) (z-W'(z)) =0 (5.7) 

Generally equation (|5.7|) has nfi solutions - the points in complex plane C, or in complexified situation (z, z) — > 
(z, z) € C 2 in the two-dimensional complex euclidean space. Equation l|5.7|l . like in the one-matrix case, can 
be then treated as a degenerate classical curve. When the Coulomb interaction between the eigenvalues is 
switched on (in (|5.4|l this is literally a two-dimensional Coulomb interaction, justifying the terminology we have 
formally introduced even in the one-matrix model case) , this set of points turns into a set of one-dimensional 
"trajectories", or into a real section of some complex curve. In what follows, we are going first to study 
the structure of this curve or Riemann surface in general position, corresponding to the complex problem in 
two-matrix model and return to a real problem later. 



5.2 Equation and genus of the curve in the complex problem 

For general polynomial potentials W and W of powers n + 1 and n + 1 correspondingly, the highest degree terms 
of the analytic equation of the classical curve 1)5. 7|) have the form z n z n + a,nZ n+1 + a n z n+1 with the coefficients 

fts = - ((« + l)*n+i) , a n = - ((n + l)t„+i) _1 (5.8) 

This means that the generic desingularization of the classical curve (|5.7() acquires the following form [2Hj (see 
also US]]) 

F(z, z) = ((n + l)(n + l^+i^+i)" 1 (z - W' n {z)) (z - W^z)) +...= 

= z n z A + a h z n+1 + a n z A+1 + ^ f ij z i z i 

(uj)e(N.p.)+ 

where we have already used complexified notations (z, z) £ C 2 for the co-ordinates embedding the curve (|5.9|) 
into the two-dimensional complex space, instead of mutually complex-conjugated co-ordinates (z, z) € C of a 
point in the eigenvalue plane of fig. El For the future purposes we will precise the form of algebraic equation of 
the curve <|5.9|l for the real potential (|2.4|) with complex-conjugated coefficients and equal powers n = h 

F(z,z) = z n z n + a n z n+1 +a n z n+1 + £ / y -zV = (51Q) 

i,je{N.p.) + 



(5.9) 
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Figure 7: The Newton polygon for the curve l|5.10|l . The highest degree terms in (|5.10() determine the shape of 
the polygon and the integer dots inside it count the number of independent holomorphic differentials, or genus 
of the curve. Clearly this number is equal to the area of "dual" square except for one (black) point, so that 
g — nh — 1 = n 2 — 1 . 



The properties of the curves l|5.9|l , l|5.10|l can be easily established via the Newton polygon (see fig. [7\ for the 
case n = h, for general values of n and h the (n + 1) x (n + 1) square on fig. [7| should be replaced by the 
rectangle of the size (n + 1) x (n + 1) with all other elements of the construction remaining intact). The first 
three terms in Ij5.9(l correspond to the three points on the boundary lines of the Newton polygon, while the sum 
over (N.P.) + in the last term stays for the sum over all points inside the Newton polygon (including the points 
on both axis not marked on fig. 0). 

The simplest basis for the holomorphic differentials on the curves (|5.9|) . H5.10J1 can be chosen as 

dvn =z l ~ Z 3— =-z l ~zi — , (5.11) 
J F z iV ( ' 

with the degrees i = i' — 1 and j = f — 1, where (i',f) £ N.P. are coordinates of the points strictly inside the 
Newton polygon, without the boundary points (see fig.[7J 12 . Counting the number of integer points inside the 
polygon one finds that the number of linear independent holomorphic differentials, or genus of the curve H5.9|l . 
equals to 

g = nh - 1 (5.12) 

For h = 1 the gaussian integration over the matrix $t returns us to the one-matrix model with the hyperelliptic 
curve of genus n — 1, considered in sect. for generic n and h the curves (|5.9(l . (|5.10|) are certainly not 
hyperelliptic. 

The generating differential, measuring the constant density of eigenvalues (|5.2|l . has the form 

dS = — ~zdz (5.13) 
2ivi 

12 For example, for n = 2 there are three points inside the polygon: > and i' + j' < 2, then the holomorphic differentials 

are labelled by i,j>0 and i + j < 1. 



27 



Figure 8: The boundary of the drop 7 and a cut of a multi-valued function z{z) inside the drop. On 7 one has 
an equality z = z(z) but this is, certainly, not true on the cut. 

The origin of this formula is demonstrated on fig. El looking at this figure one immediately comes to the following 
simple relations for the two-dimensional and contour integrals 



dz A dz = $> zdz = i> zdz = d> zdz (5-14) 

drop J*Y J~f Jcut 



meaning that one can write the eigenvalue fractions as the period integrals 

Si = / dz A dz = (f zdz (5.15) 

2m Ji_ th drop 2tti J a . 

The relations l|5.14|l allow to endow the complex curve (|5.1U|) (or (|5.9|) in the asymmetric case) with a meromor- 
phic generating differential (|5.13() . and, therefore to formulate the planar two- matrix model in the geometric 
way we have discussed above. They also clarify the relation between the complex and real problems; we are 
now going to discuss in detail the first one, intensively using the algebraic curve (|5.1U|I . and postpone discussion 
of the real problem till next section, where it will be formulated, following [3T], in terms of "doubling" fig. or 
the so called Schottky double 13 . 

The derivatives of generating differential (|5.13|) w.r.t. coefficients of the equation (|5.9|) can be computed in a 
standard way. Choosing z as a covariantly constant function, when taking derivatives over moduli of the curve, 
one writes for l|5.10|l , (|5.9|l 

FzSz + 5F = (5.16) 

where 8F = ^SfijZ 1 ^ is variation of the coefficients in equations H5.10JI . i|5.9J) . Then, for the generating 
differential (|5.13|) one gets 

S{zdz) = -SFy = -Y t ^^- ( 5 - 17 ) 

Expression (|5.17() contains decomposition of variation of the meromorphic differential (|5.13() over the basis of 
Abelian (meromorphic and holomorphic) differentials on the curve H5.1()(l . (|5.9(l . It is easy to check that the 
coefficients /y, corresponding to the meromorphic Abelian differentials of the second kind, can be expressed 
through the parameters of potential (|2.4I) of the two-matrix model, i.e. in terms of coefficients t and i of its 
harmonic part [2E|- This follows immediately from substitution into l|5.10|l . I|5.9|l of the asymptotic expansion 
of the following branch of the function z 



-1 

W\z) + O (z- 1 ) = fctfc** -1 + O (z" 1 ) (5.18) 

k=l 



13 As a demonstration of non-trivial relation between the complex and real problems, we present at fig. |5]the real section of the 
curve, described by particular cubic equation. Clearly, formula 15.121 for the complex problem gives g = 3, i.e. four independent 
variables 15.151 . corresponding to all extrema. Moreover, even the real section of complex curve in this example consists of four 
disjoint real contours. However, since cubic potential possesses a single minimum, only one of these contours can be located on 
physical sheet for a real problem. Nevertheless all extrema are indistinguishable from the point of view of complex problem and 
"holomorphic" data I5.1UI . 15.131 . The "unphysical" drops give also essential contribution to the full non-perturbative partition 
function, but the detailed discussion of this topic is beyond the scope of this paper. 
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Figure 9: Curve for 0.11z 2 z 2 — z 3 — z 3 + 6zz — 6.5 = by MAPLE computation. 



which gives rise the formulas l|5.8|) etc. Using Ij5.18|) and (|5.13|) one can immediately conclude that 



1 




(5.19) 



2nik 



The rest of the expansion l|5.17|l consists of the coefficients, corresponding to the linear combination of the third 
kind Abeian differential and the holomorphic differentials l|5.11|l . 

The dependence of the two-matrix model planar free energy upon the filling numbers (|5.15() and parameters 
of the potential l|2.4J) . as in the one-matrix case, can be geometrically formulated by 



The integrability of H5.20JI follows from the symmetry of the period matrix of the curve (|5.1()|l and the integra- 
bility of (|5.21|) from the Riemann bilinear relations analogous to l|4.7|l . 

5.3 The structure of the two-matrix model complex curve 

To understand better the structure of the curve l|5.9|) . I|5.10|l . consider first the cubic example. Writing equation 
H5.10fl for n = n = 3 with some arbitrary coefficients 




(5.20) 



where {Bi} are the canonical dual cycles A t o Bj — 5ij (see fig- EJ on the curve l|5.9|l . (|5.1U|) . and 




(5.21) 



F(z,z) 



z 2 z 2 + az 3 + az 3 + bz 2 z + bzz 2 + cz 2 + cz 2 + fzz + qz + qz + h = 



(5.22) 
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Figure 10: Cubic curve as a cover of z-plane. 



one has to make it consistent with the asymptotic (|5.18|l 



3 



z = W'(z) - G(z) = kt k z k - x + O (z- 1 ) (5.23) 
fe=i 

Substituting Ij5.23(l into the eq. I|5.22|l and collecting the coefficients of the terms z 6 , z 5 and z A one gets 



1 _ 2t 2 _ h 2t 2 8t? 

3t? ~3fc' ° ~ 3*3 ~ 9ht 3 + 9t 3 



6=4^, c=^f-^4- + 4| (5.24) 



together with their complex-conjugated counterparts, i.e. the coefficients at higher powers of equation (|5.22l) 
are indeed completely fixed by parameters of the potential l|2.4ll . Four lower coefficients /, q, q and h correspond 
to the bipole differential (|4.6(l and three holomorphic differentials and their values depend also on the filling 
fractions (|5.14|l . i.e. of the periods Ij5.15|) of the differential (|5.13|) 14 . The classical "expectation values" of these 
coefficients (at vanishing filling fractions) can be extracted from eq. (|5.7|l . 

Let us now think of the curve 15.22f> as of the Riemann surface of multi- valued function z(z), being then a 
three-sheet cover of the complex z-plane. On the first, physical sheet, there are no branching at z — > oo, as it 
follows from the asymptotic (|5.18|) . I|5.23|l . However, the "complex-conjugated" asymptotic 

z = W\z) + O (z- 1 ) (5.25) 

on two unphysical sheets z cx y/z stays, that their two infinities are glued, being an end-point of a cut, see 

fig.cni 

The branch points at z-plane are determined by zeroes of the differential dz, or by Fj = 0. Considering the 
simplest non-degenerate case of the curve (|5.22|) 

z 1 ~z 1 + az 3 + az 3 + h = (5.26) 

it is easy to see that there are nine branch points in the z-plane without infinity z — oo (of course, one comes 
to the same conclusion looking at the Cardano formula, or from the index theorem, see below). 

The structure of the curve is depicted at fig. The curve can be seen as two copies of P 1 , glued by four 
cuts, i.e. in general position it has genus (|5.12|l g = 3. There are two "infinities" z = oo, z = oo, one of 
them being a branch point. We have shown schematically on fig. 1101 the possible cuts, and the corresponding 
choice of the canonical A-cycles. In classical situation (|5.7fl the degenerate curve can be seen as two parabolas 
intersecting at four points, and under the "quantum resolution" of singularities these points turn into four cuts 
connecting two spheres P 1 as at fig. ITU1 

14 Note, that the equations of the curve 15.101 . 15.221 is often written implying some reality condition onto the coefficients, but 
as usual, the deformations of these coefficients even in this case should be considered as independent complex variables. 
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Figure 11: Generic curve of the two matrix model with symmetric potential as a cover of the z-planc. In 
contrast to fig. 1101 each fat line consists of a stack of n cuts. 

Now it is now almost clear, how the curves (|5.9|l and (|5.1U|) look for generic polynomial potentials. Say, the 
curve H5.10|l of degree n, i.e. when W'(z) ~ z n + . . . (see fig. ITT|) . can be again presented as two spheres P 1 
glued by n stacks of cuts. One of these P x 's corresponds to the "physical" z-sheet, the other one is glued at the 
oo_ from n copies of "unphysical" z-sheets. Each stack consists of n cuts, so their total number is n 2 among 
which one can choose n 2 — 1 independent, in the sense of surrounding them cycles, whose number is equal to 
the genus ()5.12(l of this Riemann surface. 

The differential dz has always a pole of the second order at oo + on the upper, or "physical" sheet, and a pole 
of the order n + 1 at oo_ since z\ oa _ oc z" + . . .. It gives altogether n + 3 poles, and from the Riemann-Roch 
theorem one concludes that the number of branching points, or zeroes of dz is equal to 

#(dz = 0) = n + 3 + 2(n 2 - 1) - 2 = 2n 2 + n - 1 (5.27) 

reproducing nine for n = 2. In general position this gives exactly 2n 2 branch points, being the end-points of n 2 
simple cuts, and n — 1 ramification points, connected by cuts with oo_. 

5.4 Degenerations of the maximal genus curve 

Up to now we have considered the complex curve of two matrix model of its maximal possible genus l|5.12|l . 
corresponding to the situation when all extrema are filled in by nonvanishing fractions of eigenvalues l|5.14|l . 
However, in many cases, as for example for the real problems to be discussed below, some part of extrema, 
corresponding to unstable configurations, can remain empty, what from the point of view of the curves l|5.9|l . 
(|5.10|l corresponds to their degenerations. Before passing to detailed discussion of how the real problem can 
be formulated in terms of the Schottky double, let us consider, following [2Hj, how the curve l|5.10[l can be in 
principle degenerated. 

The genus g — n 2 — 1 of the curve l|5.10|l decreases if there exists nontrivial solution to the following system 
of equations 

F(z, S) = 0, dF = d z Fdz + d~ z Fdz = (5.28) 

This system imposes constraints to the coefficients of of the equation (|5.1()(l . which can be found, computing 
the resultant of the equations H5.28J) . or the discriminant of the curve. However, these constraints cannot be 
effectively resolved in general position. 

To get an idea how the curve (|5.1()|l can be degenerated consider, first, the cubic case H5.22J) and let, in 
addition, all coefficients of this equation be real. Then, it is easy to see that equation ()5.22(l can be rewritten 
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Figure 12: The curve 



as double cover of the torus. When the torus 



3 curve (|5.22|l falls down to g ro d 



1. 



degenerates, the genus 



in the form 



where 15 



Y 2 + aX A + cX 2 + qX + h- - ((3a - b)X + 2c - ff = Y 2 + P(X) = 



X = z + z, 



Y = zz — - ((3a - b)X + 2c - f) 



(5.29) 
(5.30) 



The formulas l|5.30|> show that our curve l|5.10[l can be presented as a double cover of the torus (|5.29() with four 
branch points (where the transformation H5.30J1 becomes singular) being solutions to equation l|5.22|l under the 
substitution z = z. Hence, the curve (|5.22J1 can be also presented (in addition to the picture of fig. llOjl as two 
tori glued by two cuts (see fig.EJ- 

Now it becomes clear, how this picture can be degenerated. Rewriting equations l|5.28|l as 



F 2 = zF Y + P'(X) 
F- z = zF Y + P'{X) 



(5.31) 



one immediately finds that they lead either to z = z or to Fy — and, hence to P'(X) = 0. In the second 
case the torus (I5.29f) degenerates, while z = z leads to degeneration of the cover of this torus. When the torus 
degenerates into a rational curve, one gets the degeneration of the Riemann surface l|5.22|) presented as a double 
cover of sphere with two cuts, i.e. as an elliptic curve of genus g = 1 with an extra two pairs of the singular 
points (see fig. IT^ . 

Now, in the general case ()5.1(Jf) with real coefficients, the substitution analogous to (|5.3U|I brings it to the 
form 



yn , X i 







(5.32) 



where by dots we denoted monomials of lower powers in X and Y , and there are no "mixed" terms 16 . The 
genus of the curve l|5.32|l can be again easily computed by the Newton polygon (see fig. I13JI . which gives 



n(n — 1) 



(5.33) 



15 One may also "tune" for simplicity the coefficients of the potential 15.241 to get 3a = b and 2c = /. 

16 The curves of this type will be considered in the second part of this paper in the context of matrix model solutions to minimal 
string theory, or the (p, q) critical points of two-dimensional gravity with \p — q\ = 1. 
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In the same way one may present the generic curve of the two matrix model (|5.10() as a double cover of the 
Riemann surface (|5.32[l with 2n branch points. Indeed, the Riemann-Hurwitz formula 

2 - 2. 9 = # S • (2 - 2g ) - # B.P. (5.34) 

where # S is number of sheets of the cover and # B.P. is number of branch points, gives for g = n 2 — 1 and 
ffo = 9* exactly # B.P. = 2n. It means that the generic curve of the two matrix model H5.1(Jfl can be presented 
as a double cover of the curve (|5.32() with n cuts, and when the curve (|5.32l) degenerates into a rational one, 
the genus of the curve (|5.1U|I falls down to 

Sred = n - 1 (5.35) 

growing already linearly with the highest power of potential, like in the case of one-matrix model, see fig. PHI 
Finally in this section, let us say few words about the rational degenerations of (|5.10() . i.e. when its smooth 
genus completely vanishes. A particular example of such total degeneration is given by the "classical" curve 
<|5.7[1 . but the rational case can be easily studied for the generic values of coefficients in <|5.1Ull . i.e. without any 
reality restriction. 

In such situation equation of the curve l|5.10|l can be resolved via the generalized rational conformal map 

n n 

z = ™ + V^%, z=-+y^u k w k (5.36) 

k=0 k=0 

with uniformizing parameter w, and the substitution of l|5.36|l into (|5.1U|) gives a system of equations, expressing 
all coefficients in terms of parameters of the conformal map l|5.36|l . and their explicit form can be found in 
Appendix [XJ 



6 Two-matrix model and the Dirichlet problem 

In this section we continue to study the partition function of the two-matrix model (|2.5|l . It turns out, that for 
the real problem, when the maximal genus (|5.12() of the smooth curve cannot be achieved by filling only the 
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n(n-l)/2 



Figure 14: The general curve 1)5.10)1 as a double cover of the curve (j5"3l>|) with a genus 5* = nljl 2 1 - > . Similarly 
to fig. El when the curve 1)5.32)1 completely degenerates into a rational curve, the Riemann surface of the 
two-matrix model (|5.1U|) degenerates into the curve of genus g rc d = n — 1 . 




Figure 15: The deformation of domain with normal displacement n(£). 

minima of the potential (|2.4() , the most effective way to rewrite the free energy as a quasiclassical tau- function 
is based on studying the Dirichlet problem for the boundary of the eigenvalue domain, shown at fig. For the 
simply-connected domains this problem is solved in terms of the conformal maps 1)5. 36)1 . but for many eigenvalue 
drops or the multiply-connected domains the solution to the Dirichlet problem requires the whole machinery 
introduced in sect.0] 

6.1 The Dirichlet boundary problem and integrability 

Let us, first, reformulate some results of 30, 31 in the form closely connected with the quasiclassical solution 
of the two-matrix model. As we already mentioned above, the eigenvalue distribution for the real problem in 
two-matrix model 1)2.511 acquires the form of some drops in the complex plane (see fig. El, with the constant 
density 1)5. 2J) for the potential 1)2.4)1 . The density remains constant and the quasiclassical hierarchy describes now 
the deformation of the shape of these drops, which can be encoded in the solution to corresponding Dirichlet 
boundary problem. 

We start with a single drop solution, corresponding to a connected domain D, bounded by a simple smooth 
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curve. Following PUH^, consider the exterior Dirichlet problem in D c = C \ D 



u{z) = ~(f u (Od n G(z,Om (6.1) 

2lT JdD 

where G{z, z') is the Dirichlet Green function, i.e. satisfying the boundary conditions G{z, z)\ zedD = G(z, z)\ z , edD 
0, and assume that D contains the point z — 0, then D c contains the point where z — oo. Our main tool in 
this section is the Hadamard variational formula, expressing variation of the Dirichlet Green function G{z, z') 
under small deformations of domain D in terms of the Green function itself: 

SG(z,z') = ^-I d n G(z,0d n G(z',Z)6n(Z)\dZ\. (6.2) 

Here 6n(£) is the normal displacement at the boundary point £ (see fig. I15|) . and d n in l|tj.l|) and l|tj.2|) is the 
corresponding normal derivative. 

For a one-drop solution, the solution to the Dirichlet problem is equivalent [JjJ to finding a conformal map 
from D c onto the complement to unit disk \w\ > 1, or any other reference domain, where the Green function is 
known explicitly. Such conformal map w(z) exists due to the Riemann mapping theorem, then 



G(z, z') = log 



w(z) — w(z') 



w{z)w{z') — 1 



(6.3) 



where bar means complex conjugation. An example of such map is given by inverse rational function (|5.36l) . and 
this basically solves the problem for a single-drop solution of matrix model 12.5fl with a polynomial potential. 

Let tk the be moments of the domain D c — C \ D defined w.r.t. harmonic functions {z~ k /k}: 

t k =-^-[ z- k d 2 z, fc = l,2,... (6.4) 



and {ifc} be the complex conjugate moments, i.e. i k = — ^ k -f Dc d 2 zz~' k ; the coincidence of notations with 
the harmonic parameters of the matrix model potential l|2.4[) is certainly not accidental. The Stokes formula 
represents them as contour integrals, e.g. 



t k = </> z zdz = <p z zdz (6.5) 

2nik J aD 2mkJ K ' 

providing, in particular, a regularization of possible divergencies in B6.4JI and directly relating them with the 
matrix model times Ij5.19ll . if z(z) is analytic continuation of the function z(z) from dD, see fig. El Besides, we 
denote by t the area of domain D: 

1 r *. 



t = - d z z (6.6) 
I" Jd 

or the total number of eigenvalues distributed in the domain with the constant density l|5.2|) . 
The harmonic moments of D c (|6.4|) are coefficients of the Taylor expansion of the potential 

§(z,z) = -- { \og\z- z'\d 2 z' (6.7) 

7T Jo 

induced by the domain D, and directly related to the potential (|2.4(l of the two-matrix model 
V(z, z) = $(0) - $(z, z) = zz - W(z) - W(z) = \z\ 2 - J2 + ik 



,z k ) 



k>l 



(6.8) 



35 



The derivative of the potential, or resolvent (|5.6|) 



1 f d 2 z' 1 

d,$(*) = ~ / 7 = ~-G{z) (6.9) 

7T JD Z - Z n 

is continuous across the boundary and holomorphic for z £ D c while for z £ D the function + z is 
holomorphic. If the boundary is an analytic real curve, both these functions can be analytically continued, and 
there indeed exists a function z(z), analytic at least in some strip-like neighborhood of the boundary contour 
and being z(z) =zon the contour. In other words, the analytic continuation z{z) of the function z away from 
the boundary contour, directly related to the generating differential (|5.13fl . completely determines the shape of 
the boundary and is called, in this context, the Schwarz function. 

The basic fact of the theory of deformations of closed smooth curves is that the complex moments i[6.4jl . 
I|6.5jl . supplemented by real variable (|6.6[1 form a set of good local coordinates in the "moduli space" of such 
curves. Moreover, it is not an overcomplete set, as follows from the explicit construction of corresponding vector 
fields, see |3J and references therein for details. The proof of this statement is based on the observation, that 
the difference of the boundary values dtC ± (()d( of the derivative of the Cauchy integral 

C(z)dz = ^I (6.10) 
is purely imaginary differential on the boundary of D: 

9 t C(z)dz = £ I + (Cj^_Ui) da (6 . n) 

2tt« JdD V C - z (C - z) 2 J 2m J dD V (-z ) 

where C(tr; t) is paramctrization of the real curve dD(t). Hence, 

(d t c + (0 - d t c-(0) rfc = d t Cd<: - d t cdc = 2zim (6.12) 

is indeed purely imaginary, and if a i-deformation preserves all the moments tk, k > 0, the differential 
dt(d( — dtCd( extends to a holomorphic differential in D c . Indeed, if \z\ < |C| for all £ £ dD, one can expand: 

8 ( d? 00 r \ 00 

d t c + { Z )dz = ^ [ — J2 z " f D ^~ k ~ 1 ^) =Y, k ( d ^ zk ~ ldz = ° ( 6 - 13 ) 

and, since C + is analytic in D, we conclude that dtC + = 0. The expression dt(d( — dtCd( is the boundary value 
of the differential —dtC~ (z)dz which has at most simple pole at the infinity and holomorphic everywhere else 
in D c . The equality 

dtto = ^-.i (dt(d( - d t CdC) = (6.14) 

then implies that the residue at z = 00 vanishes, therefore dtC~ ' (z)dz is holomorphic. Further, any holomorphic 
differential, purely imaginary along the boundary of a simply-connected domain, must vanish in this domain 
due to the Schwarz symmetry principle on the Schottky double, obtained by attaching to D c its complex- 
conjugated copy along the boundary (from the Schwarz symmetry principle, dtC~dz extends to a globally 
defined holomorphic differential on sphere P 1 , i.e. must vanish due to the Riemann-Roch theorem). 

For a fixed point z £ D c one can consider a special infinitesimal deformation of the domain such that the 
normal displacement of the boundary is proportional to the gradient of the Green function G(z,£) (see fig. lToT) 

6 z n(0 = -ld n G(z,0 (6.15) 
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Figure 16: The elementary deformation with the base point z. 




to be called, following [2], the elementary deformations with the base point z. It is easy to see that 

5 z t = -f 6n(0\<%\ = ~ <fd n G(z,0\<%\ = e 

7TK Jqd J K 

as a direct consequence of the Dirichlet formula (|6.1f) 17 . Consider now the variation S Z X of any functional 
X = X(t) under the elementary deformation with the base point z in the leading order in e, i.e. 

dX 

^ z X = Y J ^f^tk = ^{z)X ( 6 .17) 

k 

where the differential operator V(z) is defined as 

V(z) = 3 f0+ W^ tfe + ^% fc ) (6 . 18) 



fe>i 



Note that the elementary deformation can be intuitively understood [213] as a "bump" on the boundary, continued 
harmonically into D c , see fig. El Indeed, these two are related by 5 Z oc <f \d£\d n G(z, £)<5 bump (£), and the "bump" 
deformation should be understood as a (carefully taken) limit of S z when the point z tends to the boundary 
dD. 



17 Note that the elementary deformation with the base point at oo keeps all moments except to fixed. Therefore, the deformation 
which changes only to is given by <5n(£) = — ^d n G(oo, £). 
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Fix now three points Z\, z 2) z 3 6 C \ D and compute 5 Zi G(zj, Zk) by means of the Hadamard formula 16. 211 . 
Using (|6.17() . one can identify the result with the action of the vector field V(zi) onto the Green function: 

V(z 3 )G( Zl ,z 2 ) = -±- <f d n G(zi,0d n G(z 2 ,0d n G(z3,0m (6.19) 

Remarkably, the r.h.s. of i|6.19|) is symmetric in all three arguments, i.e. 

V(zi)G(z 2 , z 3 ) = V(z 2 )G(z 3 , z x ) = V{z 3 )G( Zl ,z 2 ) (6.20) 

This is the key relation which allows to represent the deformation of the Dirichlet problem l|6.1|l as an integrable 
hierarchy of non- linear differential equations |3U) . with (|6.2U|I being the integrability condition of the hierarchy. 
It follows from (|6.2U|I that there exists a function Td — ^T>(t) of the moments Q6.4fl . (|6.6|l such that 18 

+ ^V(z)V(z')T D (6.21) 

i.e. the Green function is a double V(z)-derivative, up to a time-independent part, which is determined from 
the boundary conditions at the coinciding points z = z' and when z, z' — » oo. 

Eq. (|6.21() allows one to obtain a representation of the function Td = To(t) as a double integral over the 
domain D. Set $(2) = V(z)Td, this function is determined by its variation under the elementary deformation 

S c $(z) = -2elog|C 1 ~z- x \ +2eG((,z) (6.22) 

which is read from eq. (|6.21|> by virtue of l|6.17|l . This allows one to identify $ with the "modified potential of 
domain $(z) = $(z) — $(0) + to log \z\ 2 , where $ is given by l|6.7|l . Thus, one can write 

V(*).F D = $(*) = -- / log \z~ l - CV 2 C = «o + 2Re £ ^z" fc (6 . 23 ) 

The last equality should be understood as the Taylor expansion around infinity, and the coefficients Vk are the 
moments of domain D (the dual to l|6.4|l harmonic moments), defined as 

v h = ~ [ z k d 2 z = 8 tk T D (k > 0) , v = -$(0) = - / log \z\d 2 z = d ta T D (6.24) 
7T Jd k Jd 

i.e. the moments of the complementary domain D are completely determined by the function Td of harmonic 
moments of D c . Formulas (|6.24|l . after rewriting them as contour integrals, allow to identify the function Td = T 
with the quasiclassical tau- function corresponding to the one-drop solution of the planar matrix model i|2.5|) ■ 

In a similar manner, one arrives at the integral representation of the tau-function itself. Using the fact, that 
the elementary deformation or the operator V(£) applied at the boundary point £ G <9D (where G(z, £) = 0) 
is the bump deformation (see fig. I17fl or attaching a "small piece" to the integral over the domain D, and 
interpreting l|6.23|l as a variation S z Td we arrive at the following double- integral representation of the tau- 
function 

To = [ $(z)d 2 z = -\ f f \og\z~ 1 -C 1 \d 2 zd 2 C (6.25) 
2tt J d tt JdJd 

This is nothing, but a continuous version of the effective potential l|5.4|l . calculated on its extremal value, and 
this formula remains intact in the multi-support (or multiply-connected) case below, though an identification 
with the free energy for the multi-support solution of the two-matrix model Ij2.5|l requires slightly more care. 

18 Such formula was first conjectured by L. Takhtajan, see l.'Vil for and discussion. 



G(z, z) = log 



1 1 
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D C = C\D 




Figure 18: A multiply-connected domain D c = C \ D for g = 3. The domain D = Uq=o ^ a consists of g + 1 = 4 
disconnected parts D Q with the boundaries dD a . To define the complete set of harmonic moments, we also need 
the auxiliary points z a € D a which should be always located inside the corresponding domains. 



6.2 The multiply-connected case and the Schottky double 

Now we can turn to generic multi-support solution. Let D Q , a = 0, 1, . . . , g, be a collection of g + 1 non- 
intersecting bounded connected domains in the complex plane with smooth boundaries dD a so that D = 
U^ = qDq,, and the complement D c = C \ D becomes a multiply-connected domain in the complex plane (see 
fig. IT%|l . Let B a be the homology classes of the boundary curves <9D Q , assumed to be positively oriented as 
boundaries of D c , so that U^, =0 B a ~ — dD, or each B a has the clockwise orientation. 

Comparing to the simply-connected case, nothing is changed in posing the standard Dirichlet problem. The 
definition of the Green function and the formula 16.1|) for the solution of the Dirichlet problem through the 
Green function remain intact. A difference is, however, in the set of harmonic functions: any harmonic function 
is still the real part of an analytic function but in the multiply-connected case these analytic functions are not 
necessarily single- valued (only their real parts have to be single- valued). In other words, the harmonic functions 
may have non-zero "periods" over non-trivial cycles 19 - here over the boundary contours or B a . In general, the 
Green function has non-zero "periods" over all boundary contours, hence in the multiply-connected case it is 
natural to introduce new objects, related to the periods of the Green function. 

First, the harmonic measure w a {z) of the boundary component dD a is the harmonic function in D c such 
that it is equal to unity on 8D a and vanishes on the other boundary curves. In other words, the harmonic 
measure solves the particular Dirichlet problem, which from the formula l|6.1|l looks as 

™<x{z) = ^-l d n G(z,0\d(\, a = l,...,g (6.26) 

so the harmonic measure is the period of the Green function w.r.t. one of its arguments. From the maximum 
principle for harmonic functions it follows that < w a (z) < 1 in internal points, and moreover, it is obvious, 
that X)a=o w a( z ) = 1. In what follows we consider the linear independent functions w a {z) with a = 1, . . . , g. 

19 By "periods" of a harmonic function / we mean the integrals $ dnfdl over the non-trivial cycles. 
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Further, taking "periods" , we define 



hiD . 



Ta P = 2 f~- d n^(C)\dC\, a,p = l,...,g (6.27) 



which is a symmetric, non-degenerate and positively-definite (imaginary) matrix. It will be clear below that 
the matrix (|6.27[1 can be identified with the matrix of periods ()4.2[l of holomorphic differentials on the Schottky 
double of the domain D c . 

For the harmonic measure and the period matrix there are variational formulas similar to the Hadamard 
formula l|6.2[l . They can be derived either by a direct variation of (|6.26|) and (|6.27|) using the Hadamard formula 
or by a "pictorial" argument (see fig. 1161 ITTjl . and the formulas themselves look as 



5w a {z) = ^-i d n G(z, d n *7 a (Q 8n(Q \d£\ (6.28) 

27T JdD 



5T aP = \l d n w a (0 d n w (O 5n(0 R| (6.29) 

z JdD 



With a planar multiply-connected domain one can associate its Schottky double - a compact Riemann surface 
£, endowed with antiholomorpic involution, the boundary of the initial domain being the set of the fixed points 
of the involution. The Schottky double of the domain D c can be thought of as two copies of D c ( "upper" and 
"lower" sheets of the double) glued along their boundaries <9D C , with two infinities added (oo and 6b). In this 
set-up the holomorphic coordinate on the upper sheet is z inherited from D c , while the holomorphic coordinate 20 
on the other sheet is z. The Schottky double of the multiply-connected domain D c is a Riemann surface £ of 
genus g = #{D Q } — 1. A meromorphic function on the double is a pair of meromorphic functions /, / on D c 
such that f(z) = f(z) on the boundary <9D C , similarly, a meromorphic differential on the double is a pair of 
meromorphic differentials f(z)dz and f{z)dz such that f(z)dz = f(z)dz along the boundary curves. 

We fix on the double a canonical basis of cycles: the i^-cycles to be homologically equivalent to the 
boundaries of the holes — dD a for a = 1,. .. ,g, and the ^4 Q -cycle connects a-th hole with the distinguished 
0-th one. To be more precise, when computing integrals, fix points £ Q on the boundaries, then the A a -cyc\e 
starts from £ , goes to £ a on the "upper" (holomorphic) sheet of the double and goes back the same way on the 
"lower" sheet, where the holomorphic coordinate is z, see fig. El Harmonic w a can be represented as real parts 
vj a {z) — Lu a (z) +uj a (z) of the holomorphic multivalued functions in D c . The differentials du a are single-valued 
holomorphic in D c and purely imaginary on all boundary contours, so they can be extended holomorphically to 
the lower sheet as —duj a (z). In fact this is the canonically normalized basis of holomorphic differentials (|4.1() 
on the double S with chosen canonic basis of the cycles, since 

duip = ( dujp{z) + { (-dup(z)\ = 2Re j duj fj (z) = vap{^ a )-wp(^o) = 5 a p 
Ac Ho J£a J£o (6.30) 

i ' 



dwp = - f d n w fj (()\d(\ = T al3 

B a Z JdD a 

where the last equality gives the period matrix of the Schottky double l|6.27|l . 

One can still use harmonic moments to characterize the shape of a multiply-connected domain. However, 
the set of harmonic functions should be now extended by adding functions with poles in any hole (not only 
in Do as before) together with the functions, whose holomorphic parts are not single- valued. To specify this 
set, let us mark points z a 6 D Q , one in each hole (see fig. 1181 without loss of generality, it is convenient to 
put Zq = 0), and consider single-valued analytic functions in D c of the form (z — z a )~ k and harmonic functions 

20 More precisely, the proper coordinates should be l/z and l/z, which have first order zeros instead of poles at z = 00 (and 
z = 00). 
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So 



A V 



D c =C\D 




Figure 19: The domain D c with the A a -cycle, going one way along the "upper sheet" and back along the "lower 
sheet" of the Schottky double of D c . For such choice one clearly gets the intersection form A a o Bp = 6 a p for 
a,{3=l,...,g. 



log |l — ^| with the multi- valued holomorphic part. The arguments almost identical to the ones used in the 
simply-connected case show |31j . that the parameters 



-[ {z-z a ) n d 2 z, a = 0,1, •■■,5, n>l 

7T J D = 



together with their complex conjugate, 



- / log 




2 

d 2 z 


7T JD= 


Z 





d 2 z = $(0)-<i>(z a )-\z a \ 2 , a=l,...,g 



(6.31) 



(6.32) 



and ta oc Area(D) still given by (|6.6(l . uniquely define D c , i.e. any deformation preserving these parameters 
is trivial. A crucial step is the change of variables from M„ jQ to the variables Tk [2], which are finite linear 
combinations of the M n a 's or the moments 



t = t , T k = —<p a k (z)zdz = --l d z za k (z), k > 

2tT« JdD T , 



(6.33) 



with respect to a kind of Krichever-Novikov basis of functions (the indices a and (5 are understood modulo 
.9 + 1): 

g 



u-l 

a m ( g +i)+ a (z) = n ( z ~ z ^~ rn n ( z ~ z ^ 

0=0 0=0 
-k 



(6.34) 



3k(z) 



z 



0{z 



-k-l\ 



Any analytic function in D c vanishing at infinity can be represented as a linear combination of which is 
convergent in domains such that 1 11/3=0 — z p)\ > cons t 21 - 

21 In the case of a single hole g = formulas 16.341 give the basis used in the previous section a^ = z~ k ; note also, that always 
ao = 1, ai = l/z, therefore tq = to and ri = M i o = 
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The parameters r k , 4> a can be treated as local coordinates in the space of multiply-connected domains, 
analogously to the simply-connected case, the details of the proof can be found in |31|. see also Appendix iBl 
Instead of <p a one can use the already known fractions of eigenvalues, or the areas of the holes 

Area(D Q ) 1 f , 2 1 / _ , 1 / _ , „ 

S a = = — a z = <t> zdz = <t> zdz , a = 1, ... ,g (6.35) 

T T JD a 2m JdD a 2lTl JBo, 

where we have again introduced the Schwarz function z{z) = z in some strip-like neighborhoods of all 

zedD 

boundaries. Note, that in the notations of this section the ^-variables are expressed naturally through the 
_B-pcriods of the generating differential, in contrast to the more common choice l|5.15|l used before. This is 
related to the fact, for the Schottky double it is convenient to choose the canonical differentials l|6.3U[l related 
to harmonic measure and normalized to the ^-cycles, as their are chosen on fig. 1191 For these different choices 
the tau-functions are related by the duality transformation, already mentioned above and to be discussed in 
detail in sect. 16.51 

At the same time, the variables <p a do not have any geometric sense on the Schottky double, instead, one 
can introduce the variables U a (again the coincidence of notations with the previously introduced Lagrangian 
multipliers is not accidental), which were referred to in as virtual A-periods of the differential z(z)dz on 
the Schottky double, i.e. coinciding with these yl-periods always, when the latter can be rigorously defined. 
Following |3J , consider the basis of differentials db^ , satisfying the "orthonormality" relations 22 

— / a k db k , = 5k,k' (6.36) 
Jki Jq D 

for all integer k,k' 6 Z. Explicitly they are given by: 

dz 9 

db m(g+1)+a = —— J[(z- z ) m J[(z- zp-x) (6.37) 

9 0=0 0=0 

where we identify Z-\ = z g . Now one can introduce the meromorphic differential on £ with the only pole at oo 
on the upper sheet, where it has the form 

dh k {z) = db k (z) + 0{z~ 2 )dz (6.38) 

and vanishing yl-periods 

j> d£l k =0 (6.39) 

i.e. it is a canonically normalized meromorphic differential with a singular part H6.38(l . The normal displacements 
of the boundary given by real and imaginary parts of the normal derivative d n £l k define a complex tangent vector 
field (the partial derivatives at constant IT a ) 

9 "=^-5>(z a )^ (6>40) 

Q 

to the space of multiply-connected domains, where 

b k {z)= f db k (6.41) 
Jo 

is a polynomial of degree k. These vector fields keep fixed the formal variable 

n Q = 4> a + 2Re^ b k (z a )T k (6.42) 



2 Some properties of the functions 16.341 and their dual differentials 16.371 are collected in Appendix IBl 
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The sum in the r.h.s. generally does not converge, but in the case when the Schwarz function has a meromorphic 
extension to the double E, the sum converges and the corresponding quantity coincides with the ^-period of 
the extension of the differential zdz, see formula l|tj.71fl below. 

As in the simply-connected case, we introduce the elementary deformations 

6 Z with Sn(0 = -^d n G(z,0, z E D c 

e ( 6 - 43 ) 
<5 (q) with tfn(£) = - -d n w a (£) , a=l,...,g 

where w a (z) is the harmonic measure of the boundary component <9Dq, (see <|6.26(l '). Then the variations of the 
local coordinates under elementary deformations are: 



S z T k = ea k (z), 5 z <j> a = elog S^r k = 0, S^4>p = -2e5 a p (6.44) 



since 



S, [ /(C) A = 5 / f(()d n G(z,0\dC\ = -enf(z) 

JD C 1 JdD c 



(6.45) 



for any harmonic function / in D c , and 



6 {a) / f{Qd 2 t = -& fd n vj a \dC\=-- & d n f\dC\ = -ie& d c fd( (6.46) 

leading to (|6.44|l for particular choices of /. Variations of the variables II a (in the case when they are well- 
defined) then read 

5 z Tl a = 0, 6^U = ~2eS aP (6.47) 
and for any functional X on the space of the multiply-connected domains the following equations hold 

6 Z X = eV(z)X 

, n (6.48) 

S (a) x = _ 2e d*X = -2ed%X V ; 

The differential operator V(z) is defined now by the formula 

v(z) = dl + £ (* k {z)dl +i^I)a T n fc ) (649) 

k>l 

The functional X can be regarded as a function X = X^(4> a ,T k ) on the space with the local coordinates 4> a , T k , 
or as a function X = X n (H a , r k ) on the space with the local coordinates n Q , r k (we stress again, that although 
the variables II Q are generally formal, their variations under elementary deformations and the vector-fields d^ k 
are always well-defined). 

Let K(z,()d( be a unique meromorphic Abelian differential of the third kind 14. 5|) . I|4.6|l on the Schottky 
double E with the simple poles at z and oo on the upper sheet with residues ±1. Then 

2d c G{z, ()d( - 2<9 c G(oo, ()d( = K{z, ()d( + K(z, ()d( (6.50) 
and the differential d(l k (() H6-38fl . (|6.39|) can be represented in the form 

dn k (C) = ^-j K ^ Odb k (u) (6.51) 

where the w-integration goes along a big circle around infinity. Using summation formulas (|B.1|) (see Appendix |B|) 
we obtain that 

- £ B k (z)dQ k (0 = <f K(uX)du = R ^ ^ 

fc>l Lm J°o U Z 
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Therefore, for the Green function one gets an expansion 

2d c G(z, ()d( = cM (0 - (a fc (aOdfi fc (0 + c.c.) (6 53) 

fc>i 

where the complex conjugated d£Lk is a unique meromorphic differential on £ with the only pole at infinity 6b 
on the lower sheet with the principal part — dbk{z) and vanishing A-periods. 

6.3 Quasiclassical tau-function of the multi-support domain 

Applying the formulas (|6.2|l . i|6.28[l . (|6.29|) . one finds the variations of the Green function, harmonic measure 
and period matrix under the elementary deformations: 

5 Zl G(z2, z 3 ) = S Z2 G(z 3 , zi) = S Za G(zi,z 2 ) 
5 Zl w a (z 2 ) = 8 ( - a) G(z 1 ,z 2 ) = 8 Z2 w a (zi) 

6l oi >wfi{z)=6Wm a {z) (6.54) 
S z T a0 = i^wpiz) 

From 16.48f) it follows that formulas l|6.54[) can be rewritten in terms of the differential operators V(z) and 
d a = d/d(j) a = d/dU a : 

V( Zl )G(z 2 ,z 3 ) = V(z 2 )G(z 3 , z x ) = V(z 3 )G(z!, z 2 ) 
V{z 1 )w a {z 2 ) = -2d a G(z 1 ,z 2 ) 

d a vjp(z) = dpm a {z) (6.55) 
V(z)T Q/ 3 = — 2irid a 'Ui7p(z) 
daTpj = dpT ia = d^Tap 

These integrability relations generalize formulas (|6.2()(l to the multiply-connected case; the first line just coincides 
with H6.20J) . while the other ones extend the symmetricity of the derivatives to the harmonic measure and the 
period matrix. 

Again, (|6.55|l can be regarded as a set of compatibility conditions of an infinite hierarchy of differential 
equations. They imply that there exists a function To = JFd(1I q ,t) such that 

G( Zl ,z 2 ) ^log^r 1 ^- 1 ! + iv(zi)V(z 2 )^D 

m a {z) = -d a V{z)T D ( 6 - 56 ) 
T af 3 = 2ni d a dpT D 

The function To is the (logarithm of the) tau-function of multiply-connected domains and will be related below 
in sect. !6.5l with the free energy of multi-support solutions to the matrix model (|2.5|) after duality transformation. 

Set again <t(z) = V(z)To, equations l|6.56fl then determine &(z) for z e D c via its variations under the 
elementary deformations 

6 c $(z) = -2e\og\C 1 -z- 1 \+2eG((,z) 
6 {a) $(z) = 2ezu a (z) 

Indeed, using H6.43J) . for the variation of 

*{?) = -- [ \og\z- l -C 1 \d 2 C = ^)-^) + r \og\z\ 2 (6.58) 
~x Jo 
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just coinciding with (|6.23|) . if D is understood as union of all D a 's, one gets 

k{~~ I log\z- 1 -z'- 1 \d 2 A=-<f \dt\d n G(C,0\og\z- 1 -C 1 \ = 

\ K JD J K J dD c 

= - I \dt\d n G(C,0 (log Iz- 1 - r X | - G(z,0) = -2elog Ir 1 - z- 1 ] + 2eG(C, z) 
n Jao c 



Similarly, for z£ D c one obtains 



e 



\d£\d n (log\z- 1 -C 1 \-G(z,0)=2ew a (z) 

n 'dD a 



The same calculation for z € D yields 



The coefficients of an expansion of $ at infinity define the dual moments v k : 

V{z)T D = *(z) = -- / loglz- 1 - CVC = vo + 2Re V v k a k (z) 
7T In 

JD k>Q 

which arc moments of the union of the interior domains with respect to the dual basis 



(6.59) 



= -/ |de|fl^w a (0(log|a- 1 -r 1 |-G(*,0)= (6-60) 

J9D C 



I 2ed Q/3 n z £ Dp , /3=l,...,g 



(6.62) 



^fc = - I b k (z)d 2 z (6.63) 

From IpTH!^ it follows that 

v k = d? h F D (6.64) 
and the same arguments show that the derivatives 

S a = - d a T D (6.65) 

are just areas of the holes (|6.35|) . Indeed, 



(6.66) 



The integral representation of To is found similarly through its variations, read from l|6.62|) . and the result is 
given by the same formula l|6.25|) as in the simply-connected case provided by D = U^ =0 D a is now the union of 
all D a 's. 
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6.4 Polynomial potentials and algebraic domains 



The domains are called algebraic if the Schwarz function z(z) has a meromorphic extension to a Riemann 
surface £ with antiholomorphic involution, then £ can be naturally divided in two "halves" (say upper and 
lower sheets) exchanged by this involution. The domain D c is algebraic if and only if the Cauchy integrals 

C a (z) = ^-I forzeD a (6.67) 

2?" Jqd C ~ z 

are extendable to a rational meromorphic function S(z) (the same for all a!) on the whole complex plane C 
with a marked point at infinity. Equality z(z) = S(z) — C~(z), valid by definition for z £ <9D C , can be used for 
analytic extension of the Schwarz function. Since C~(z) is analytic in D c , S(z) and z{z) have the same singular 
parts at their poles in D c . One may treat z(z) as a function on the Schottky double extending it to the lower 
sheet as z. 

It is also convenient to introduce the Abelian integral 

b(z) = [ S(z)dz (6.68) 



o 



which is multi-valued if S(z) has simple poles (to fix a single-valued branch, one has to make cuts from oo to 
all simple poles of S(z)). In neighborhoods of the points z a £ D a , see (16.34(1 . 1(6.3711 . one has 

S{z)dz = ^rfcdb fc (z), b(z) = ^T fe b fe (z) , g QQ s 

k>l k>l 

Formula ((6.42|) shows that for the algebraic domains the variables IIq, , introduced in the general case as formal 
quantities, are well-defined and equal to the A-periods of the differential z(z)dz on the Schottky double £. 
Indeed, using the fact that C (z) and C a (z) represent restrictions of the same function S(z), one can use 
representation ((B.7() from Appendix iBl and rewrite it in the form 

2 Re I / S{z)dz+ [ (z + C~ (z) - S(z))dz\ (6.70) 



J to J 

Combining this equality with the definition of II Q 1(6.42(1 . we obtain: 

II Q , = 2Re/ (z(z)-z)dz= [ (z(z)dz - zdz(z)) = I zdz (6.71) 

ho ho v ' JA a 

As an example of algebraic domains, it is instructive to consider the case with only a finite number of the non- 
vanishing moments Tk, i.e. tu — for k > n+ 1, which corresponds to the two-matrix models with polynomial 
potentials ((2.4(1 . Then z (z) extends to a meromorphic function on £ with a pole of order n at oo and a simple 
pole at 6b. The function z extended to the lower sheet of the Schottky double as z(z) has a simple pole at oo 
and a pole of order n at do. For such a domain D c the moments with respect to the Laurent basis 1(6.4(1 . as in 
the simply-connected case, coincide with the coefficients of the expansion of the Schwarz function near oo: 

n+l 

z(z) = kt k z k - x + Oiz- 1 ), z^oo (6.72) 

k=l 

literally coinciding with formula ((5.18(1 . The normal displacement of the boundary of an algebraic domain, 
which changes the variable t k keeping all the other moments (and II Q ) fixed is defined by normal derivative of 
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the function 2 Re J z dSlf., where dVL k is a second kind Abelian normalized meromorphic differential (|4.4|l on E 
with the only pole at oo 

dVL k =d(z k + 0(z- 1 )), j dfl k =Q (6.73) 

Let also d£lo be a third-kind Abelian differential (|4.5|) . (|4.6|l on E with simple poles at two infinities oo and do, 
its Abelian integral 

logw(z) = / dn (6.74) 

defines in the neighborhood of oo a function w(z) which has a simple pole at infinity. The dependence of the 
inverse function z(w) on the variables t k is described by the Whitham equations for the two-dimensional Toda 
lattice hierarchy of the form 

dlz{w) = {fl k (w),z(w)} = ^M 9t0 zH - dMw)^^ (6.75) 

In this case one can also write a quasi-homogeneity condition |2.'{j for the quasiclassical tau-function, which 
acquires the form 

1 1 9 

2-F D = - ^ To + t v + 2 ~ k )( T k"k + f k D k ) - n«S a (6.76) 

fe>l a=l 

Algebraic domains of a more general form correspond to the universal Whitham hierarchy |23| . the discussion 
of this issue can be found in 



6.5 The duality transformation and free energy of two-matrix model 

Passing from II Q to S a is a particular case of duality transformation l|4.17|l which is equivalent to the interchang- 
ing of the A and B cycles on the Riemann surface E, see fig- 01 For the quasiclassical tau-function, depending 
also upon some extra r- variables this is achieved by the partial Legendre transform (n a , t) — > T{S a ,T), 
where 

9 

T D = T+^T\ a S a (6.77) 

a=l 

The function T is dual prepotential to the tau-function of the Dirichlet problem, it solves the modified Dirichlct 
problem and can be identified with the free energy of two-matrix model in the planar large N limit in the case 
when the support of eigenvalues consists of a few disconnected drops. 

The main properties of F follow from those of Tq. According to l|6.64ll . (|6.65() we have dTo = — S a dTL a + 
J2k v kdr k (for brevity, k is assumed to run over all integer values, r^ k = f k ), so dT = J2 a Tl a dS a + J2 k v kdT k . 
This gives the first order derivatives: 

The second order derivatives are transformed as follows (see e.g. set 

_ _ d 2 T D _ d 2 T D _ _ d 2 T D 
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and similarly for T , then 

= )a0 
J~Dak = / X -F" )a^F~ jk 

T =i (6.80) 
g 

7,7'=1 

where (^ r_1 ) Q /3 means the matrix element of the matrix inverse to the g x g matrix T a p. 

Using these formulas, it is easy to see that the main properties (|6.56[1 of the tau-function are translated to 
the dual tau-function as follows: 

G{z, C) = log - C 1 1 + l -V{z)V{QT 

2iri w a {z) = - d s y{z)J : (6.81) 

where Tfc-derivatives in V(z) are taken at fixed S a . The objects in the left hand sides of these relations are: 



G(z, C) = G(z, () + in ^2 Va(z)f a p <Jp(Q 



a, 0=1 

» - (6.82) 

io a (z) = } T a pup(z) 

0=1 
f = 

The function G is the Green function of the modified Dirichlet problem, see The matrix T is the matrix 
of A-periods of the holomorphic differentials doj a on the double £ (so that zb a (z) = Cb a {z) +Q a (z)), normalized 
with respect to the B-cycles § B dujp = —8 a p, § A dibp — T a p, i.e. more precisely, the change of cycles under 
duality transformation is A a — > B ai B a — > — A a . 

An important fact is that by a simple rescaling of the independent variables one is able to write the group of 
relations (|6.81|l for the function T in exactly the same form (|6.56|l . so that they differ merely by notation. There- 
fore we will not distinguish between these two cases in |34| when studying equations, satisfied by quasiclassical 
tau-functions of matrix models. 



7 Conclusion 

In the second part of this paper [31 we consider certain examples of applications of the above methods and dis- 
cuss a similar quasiclassical geometric picture, arising already in the context of multidimensional gauge theories 
and the AdS/CFT correspondence. We also planning to discuss shortly the open problems and speculations 
how this picture can be treated beyond the quasiclassical limit. 
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NWO project 047.017.015, the ANR-05-BLAN-0029-01 project "Geometry and Integrability in Mathematical 
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Appendix 



A Rational degenerations 

Substituting (|5.36|) into (|5.10() and computing the residues one finds that the expressions 

Ri[F] = res (^w l F(z{w),z(w)) \ = (A.l) 

for / = —n(n + f ), . . . , n(n + f ) form a triangular system of equations onto the coefficients /y. It means that 
each of the equations (|A. 1|) is linear in one of the coefficients, and can be resolved step by step, starting from 
the ends of the chain. 

For the cubic potential (n = 2) the solution is 

r 2 

a = 

b = 2 uo 

Uituq _ 2 _ r 2 u r 3 u x 

c = h u — 2 r u\ + a 

m 2 u 2 u 2 u 2 

2 7,4 ru\UQ r 2 U\Ui ru^Ui 

f = r — 2 U2 U2 + 4 Mo Mo + u\ u\ — 2 — = 1 = 2 

u 2 u 2 u 2 u 2 U2 u 2 

2 2 2 -2 4 - 

r uf. _ 2 „ t u\ r uo 

q = —3 h 2U2U2UQ — U2U X — 2 — — 3 r M2 ui + Mo u\ u\ 

u 2 ' ~ u 2 u 2 u 2 

., - -2 ruiwi 2 2 - , o r-uiM Uo . „ r 3 u x 

+4 Mo r ui — 2 m u o 7" Mo + 2 h 3 

U2 "2 u 2 (A.2) 

„-,2 „,2„-: „, 



U 2 U 2 U 2 

r 6 r 2 uf. , MorMiMn 

ft = + — T- 51 + U2 M Mi 2 - _ 

M 2 M 2 M 2 M 2 

U2ru\ r 2 Mo 3 r 4 Mi Mi u Q u Uir 2 u 
- -\- -|- 2 ~\- 



M 2 M 2 M 2 M 2 M 2 M 2 

„, 2 „-; „, „. 



«2 

M Mo r4 2-2 MiTMqMo 2 

-Ml Ml M M H = h M2 M 2 



"2 


u 2 




, "0 Ml 






M 2 M 2 






+ M2 M M 2 - 


r 3 m mi 


m 2 Mi 3 r 

U2 


, r 3 Mi M§ 


r 3 Mi Mo 2 


2 -2 2 
Ml Mi r z 


M 2 M 2 


U 2 U 2 


M 2 M 2 


2 2 — 
Ml r Mo 

— + 


uxujruo 

U-2 


r 2 M M 


— 2 r Mi mq 2 - 


, o w r 2 M 2 

«2 


r 3 M Mi 
"2 



M 2 M 2 ' M 2 

— Ml Ml M2 U2 — 3 r 2 M 2 M2 + 3 M2 r Ml Mo — 2 M2 M2 Mo Mo — Ml r 2 U\ + 3 Mo V U 2 U\ 

together with the "complex conjugated" expressions for a, b, c and q, where one should replace Uk by Uk and 
vice versa. Resolving (|A. II one gets the explicit description of the rational degeneration of the curve l|5.10|l in 
terms of the coefficients of conformal map ()5.36|) . However, in general situation they are only implicitly defined 
through the parameters of the potential V(z, z). The simplest example of such degeneration is given by 



z 2 z 2 - Kz 3 - —z 3 + ( t H = - mht I zz + 

M 3i V 9tt ' 



3i§(l - 3tt) - 27t 3 tt - = (A.3) 



and this equation can be resolved via conformal map 

+ z = ^+3t to 2 (A.4) 



w 2 w 
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B Co-ordinates in multiply-connected case 



The existence of a well-defined dual basis of differentials (|6.37|) . obeying the orthonormality relation Ij6.36|l . is 
the key feature of the basis functions Zk (|6.34|) . which makes their moments ()6.33(l good local coordinates, 
in contrast to -M„ jQ , since for the functions (z — z a )~ n one cannot define the dual basis. 

The summation formulas 



^ = E dCa„(C)db n (z), | ]J0- zp)\ < | I](C - 

' n=l /3=0 /3=0 

d df °° 9 9 

£-1 = - £ dCa_„(C)db_„(z), I J] (z - ^)| > I J] (( - zp)\ 

' Z n=0 /3=0 j3=0 



(B.l) 



which can be checked directly allow us to repeat arguments of sect. 16. ll Indeed, the Cauchy integral Hfi.lOjl 

C(z)d Z = £ / (B.2) 

where the integration now goes along all boundary components, defines in each of the holes D a analytic dif- 
ferentials C a (z)dz (analogs of C + (z)dz in the simply-connected case). In the complementary domain D c the 
Cauchy integral still defines the differential C~ (z)dz holomorphic everywhere in D c except for infinity where it 
has a simple pole. The difference of the boundary values of the Cauchy integral is equal to z: 

C a (z) -C-(z) = z, z G dD a (B.3) 

From equation (|6.11() . which can be written separately for each contour, it follows that the difference of the 
boundary values of the derivative of the Cauchy integral l|B.2|l 



(d t C a (0-d t C-(())dC (B.4) 

is, for all a, a purely imaginary differential on the boundary dD a . The expansion IjB.lJl of the Cauchy kernel 
implies that if a t-deformation preserves all the moments Tfc, k > 0, then dt(d( — dtC,dC, extends to a holomorphic 
differential in D c . 

Indeed, for z close enough to any of the points z a , one can expand dtC a {z) for each a as 

-.00 f „ \ 00 

d t C a (z)dz = —J2db k d t I J> a fc (C)CrfCj = Y, dtTk dbk ^ ( B - 5 ) 

1:1 k=l \J0D / k=1 

and conclude that it is vanishes identically provided d t T^ = 0. Hence d t C~ (z)dz is the desired extension of 
d t (d( — d t (d( and has no pole at infinity due to d t r = 0. 

Using the Schwarz symmetry principle we extend d t C~ (z)dz to a holomorphic differential on the Schottky 
double. If the variables S a ()6.35(l are also preserved under the t-deformation, this holomorphic differential has 
vanishing periods along all cycles B a and therefore, it identically vanishes. This completes the proof of the 
statement, that any deformation of the domain preserving all and S a is trivial. In this proof the variables 
S a were used only at the last moment, to show that extension of dtC~ (z)dz as a holomorphic differential on 
the Schottky double E is trivial; instead the variables <fi a Ij6.32|l can be used in a similar way. 

Let us show that if they are preserved under t-deformation then the A a -periods of the extension of dtC~ (z)dz 
vanish, and therefore this extension is identically zero. Indeed, the variable 4> a ()6.32|) can be represented in the 
form 

2 r Za r d 2 C 
Q = Re / dz / —±- (B.6) 

77 Jo Jd< z - C 
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The differential ^ £_£. i s equa i to C a (z)dz for z S D Q and (z + C~(z))cfe for z G D c . Let £ , £a be the 
points where the integration path from z — to z = z a intersects the boundary contours Bq, B a . Then 

cj) a = -2 Re C°(z)dz + " C a {z)dz + (z + C~(z))dz^ (B.7) 

It is already shown above that if a i-dcformation preserves the variables then all dtC a (z)dz = 0. Thus 
vanishing of the i-derivative dt<j) a = implies 

= -dt<j) a = 2 Re I d t C~(z)dz (B.8) 

The r.h.s of this equation is just the >l Q -period of the holomorphic extension of the differential dtC~ (z)dz. 
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